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Abstract. The phenomena implied by the existence of quantum vacuum fluctuations, 
^^ , grouped under the title of the Casimir effect, are reviewed, with emphasis on new 

results discovered in the past four years. The Casimir force between parallel plates is 
^ I rederived as the strong-coupling limit of J-function potential planes. The role of surface 

\l I divergences is clarified. A summary of effects relevant to measurements of the Casimir 

^^i^ ' force between real materials is given, starting from a geometrical optics derivation of 

\^^ . the Lifshitz formula, and including a rederivation of the Casimir-Polder forces. A great 

CP ] deal of attention is given to the recent controversy concerning temperature corrections 

to the Casimir force between real metal surfaces. A summary of new improvements 

to the proximity force approximation is given, followed by a synopsis of the current 

experimental situation. New results on Casimir self-stress are reported, again based 

JL ■ on (5-function potentials. Progress in understanding divergences in the self-stress of 

g_) I dielectric bodies is described, in particular the status of a continuing calculation of the 

i-Q ■ self-stress of a dielectric cylinder. Casimir effects for solitons, and the status of the 

so-called dynamical Casimir effect, are summarized. The possibilities of understanding 
dark energy, strongly constrained by both cosmological and terrestrial experiments, in 
5^ ' terms of quantum fluctuations are discussed. Throughout, the centrality of quantum 

5^ I vacuum energy in fundamental physics in emphasized. 
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1. Introduction 

The essence of quantum physics is fluctuations. That is, knowing the position of a 
particle precisely means losing all knowledge about its momentum, and vice versa, and 
generally the product of uncertainties of a generalized coordinate q and its corresponding 
momentum p is bounded below: 

AqAp>^, (1.1) 

which reflects the fundamental commutation relation 

[q,p]=th. (1.2) 

The Hamiltonian commutes with neither q nor p in general; this means that in an energy 
eigenstate the fluctuations in q and p are both nonzero: 

Ag > 0, Ap> 0. (1.3) 

Moreover, a harmonic oscillator has correspondingly a ground-state energy which is 
nonzero: 

Ei,o,n = fj^ (n + -]. (1.4) 

The apparent implication of this is that a crystal, which may be thought of, roughly, as 
a collection of atoms held in harmonic potentials, should have a large zero-point energy 
at zero temperature: 

atoms 

uj being the characteristic frequency of each potential. 

The vacuum of quantum field theory may similarly be regarded as an enormously 
large collection of harmonic oscillators, representing the fluctuations of, for quantum 
electrodynamics, the electric and magnetic fields at each point in space. (Canonically, 
the momentum-coordinate pair correspond to the electric field and the vector potential.) 
Put otherwise, the QED vacuum is a sea of virtual photons. Thus the zero-point energy 
density of the vacuum is 

u = ylnw = 2 [ ^hc\ki (1.6) 

where k is the wavevector of the photon, and the factor of 2 reflects the two polarization 
states of the photon. 

This is an enormously large quantity. If we say that the largest wavevector 
appearing in the integral is K, say UcK ~ 10^^ GeV, the Planck scale, then U ~ 10^^^ 
GeV/cm^. So it is no surprise that Dirac suggested that this zero-point energy be simply 
discarded, as some irrelevant constant P] (yet he became increasingly concerned about 
the inconsistency of doing so throughout his life j2j). Pauli recognized that this energy 
surely coupled to gravity, and it would then give rise to a large cosmological constant, 
so large that the size of the universe could not even reach the distance to the moon 
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[Slllj. This cosmological constant problem is with us to the present JHIEI- But this was 
not the most perplexing issue confronting quantum electrodynamics in the 1930s. 

Renormalization theory, that is, a consistent theory of quantum electrodynamics, 
was invented first by Schwinger j7j and then Feynman jH] in 1948; yet remarkably, 
across the Atlantic, Casimir in the same year predicted the direct macroscopically 
observable consequence of vacuum fluctuations that now bears his name jH] . This is the 
attraction between parallel uncharged conducting plates that has been so convincingly 
demonstrated by many experiments in the last few years jTO]. Lifshitz and his group 
generalized the theory to include dielectric materials in the 1950s [TTl^l^lE]- There 
were many experiments to detect the effect in the 1950s and 1960s, but most were 
inconclusive, because the forces were so small, and it was very difficult to keep various 
interfering phenomena from washing out the effect ^Hl- However, there could be very 
little doubt of the reality of the phenomenon, since it was intimately tied to the theory of 
van der Waals forces between molecules, the retarded version of which had been worked 
out by Casimir |[TS| just before he discovered (with a nudge from Bohr [T2j) the force 
between plates. Finally, in 1973, the Lifshitz theory was vindicated by an experiment 
by Sabisky and Anderson [TH] . 

But by and large field theorists were unaware of the effect until Glashow's student 
Boyer carried out a remarkable calculation of the Casimir self-energy of a perfectly 
conducting spherical shell in 1968 ^Hl- Glashow was aware of Casimir 's proposal [2IIj 
that a classical electron could be stablized by zero-point attraction, and thought the 
calculation made a suitable thesis project. Boyer's result was a surprise: The zero- 
point force was repulsive for the case of a sphere. Davies improved on the calculation 
pT] : then a decade later there were two independent reconfirmations of Boyer's result, 
one based on multiple scattering techniques P2j (now undergoing a renaissance, for 
example, see [22]) and one on Green's functions techniques [211 (dubbed source theory 
PH] ) . Applications to hadronic physics followed in the next few years [SEl 123 I2HJ , and 
in the last two decades, there has been something of an explosion of interest in the field, 
with many different calculations being carried out [221 QH] ■ 

However, fundamental understanding has been very slow in coming. Why is the 
cosmological constant neither large nor zero? Why is the Casimir force on a sphere 
repulsive, when it is attractive between two plates? And is it possible to make sense of 
Casimir force calculations between two bodies, or of the Casimir self-energy of a single 
body, in terms of supposedly better understood techniques of perturbative quantum 
field theory [^? As we will see, none of these questions yet has a definitive answer, yet 
progress has been coming. Even the temperature corrections to the Casimir effect, which 
were considered by Sauer [31], Mehra [32], and Lifshitz [TT| in the 1950s and 1960s, have 
become controversial [S3IMllS3EaiSZllSHllSnilinillIlll2lll3llllll3linillI|- Thus 
recent conferences on the Casimir effect have been quite exciting events |I^ EH] . It is 
the aim of the present review to bring the various issues into focus, and suggest paths 
toward the solutions of the difficulties. It is a mark of the vitality and even centrality 
of this field that such a review is desirable on the heels of two significant meetings on 
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the subject, and less than three years after the appearance of two major monographs 
[TUl |2ni on Casimir phenomena. There are in addition a number of earher, excellent 
reviews |SII1EI1E21, as well as more specialized treatments ISHl EH EHl EEl • Throughout 
this review Gaussian units are employed. 

This review is organized in the following manner. In section |21 we compute Casimir 
energies and pressures between parallel 6 function planes, which in the limit of large 
coupling reproduce the results for a scalar field satisfying Dirichlet boundary conditions 
on those surfaces. Although these results have been described before, clarification of the 
nature of surface energy and divergences is provided. TM modes are also discussed here 
for the first time. Then, in section El we rederive the Lifshitz formula for the Casimir 
force between parallel dielectric slabs using a multiple reflection technique. The Casimir- 
Polder forces between two atoms, and between an atom and a plate, are rederived. 
After reviewing roughness and conductivity corrections, a detailed discussion of the 
temperature controversy is given, with the conclusion that the TE zero-mode absence 
must be taken seriously, which will imply that large temperature corrections should be 
seen experimentally. New approaches to moving beyond the proximity approximation 
in computing forces between nonparallel plane surface are reviewed. A discussion of the 
remarkable progress experimentally since 1997 is provided. In section H] after a review 
of the general situation with respect to surface divergences, TE and TM forces on 6- 
function spheres are described in detail, which in the limit of strong coupling reduce 
to the corresponding flnite electromagnetic contributions. For weak coupling, Casimir 
energies are finite in second order in the coupling strength, but divergent in third order, 
a fact which has been known for several years. This mirrors the corresponding result 
for a dilute dielectric sphere, which diverges in third order in the deviation of the 
permittivity from its vacuum value. Self-stresses on cylinders are also treated, with a 
detailed discussion of the status of a new calculation for a dielectric cylinder, which 
should give a vanishing self-stress in second order in the relative permittivity. Sectional 
briefiy summarizes recent work on quantum fiuctuation phenomena in solitonic physics, 
which has provided the underlying basis for much of the interest in Casimir phenomena 
over the years. Dynamical Casimir effects, ranging from sonoluminescence through the 
Unruh effect, are the subject of section IHl The presumed basis for understanding the 
cosmological dark energy in terms of the Casimir fluctuations is treated in section 
where there may be a tight constraint emerging between terrestrial measurements of 
deviations from Newtonian gravity and the size of extra dimensions. The review ends 
with a summary of perspectives for the future of the fleld. 

2. Casimir Effect Between Parallel Plates: A (5-Potential Derivation 

In this section, we will rederive the classic Casimir result for the force between parallel 
conducting plates p. Since the usual Green's function derivation may be found in 
monographs ^\, and was recently reviewed in connection with current controversies 
over flniteness of Casimir energies ^3) ^^ ^ill here present a different approach, based 
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on (5-function potentials, which in the hniit of strong couphng reduce to the appropriate 
Dirichlet or Robin boundary conditions of a perfectly conducting surface, as appropriate 
to TE and TM modes, respectively. Such potentials were first considered by the Leipzig 
group j2Hl EH], but recently have been the focus of the program of the MIT group 
[nm EH innj ■ The discussion here is based on a recent paper by the author [021 • ^^ ^^^t 
consider two 5-function potentials in 1 + 1 dimensions. 



2.1. 1 + 1 dimensions 

We consider a massive scalar field (mass /x) interacting with two (5-function potentials, 
one at x = and one at x = a, which has an interaction Lagrange density 



C, 






(2.1) 



where we have chosen the coupling constants A and A' to be dimensionless. (But see 
the following.) In the limit as both couplings become infinite, these potentials enforce 
Dirichlet boundary conditions at the two points: 

A,A'^oo: 0(0), 0(a) ^0. (2.2) 

The Casimir energy for this situation may be computed in terms of the Green's 
function G, 



G{x,x') = i{T^{x)4>{x')) 

which has a time Fourier transform, 

duj 



G{x,x') 



27r 



e-''^^'-''^g{x,x';u;) 



(2.3) 



:2.4) 



Actually, this is a somewhat symbolic expression, for the Feynman Green's function ()2.3j] 
implies that the frequency contour of integration here must pass below the singularities 



in u on the negative real axis, and above those on the positive real axis 
reduced Green's function in (12. 4p in turn satisfies 



o 9 A ^, , A ^, 

■TTT + K^ + -6{x) + —6{x 
ax^ a a 



g{x,x') = 6{x — x'). 



The 



(2.5) 



Here k^ = ^^ — uo"^. This equation is easily solved, with the result 



g[x,x 
A 



1 ^-k\x-x'\ j^ 1 



2/t 



2k A 



A' 



2K,a —K,{x+x') 



w 

(2/ta)2 
A' 



2cosh/T;|a; — x' 
A 



1 



K,{x+x') 



2Ka \ 2Ka J 2Ka \ 2Ka 

for both fields inside, < x,x' < a, while if both field points are outside, a < x, x', 



(2.6a) 



g{x,x') = —e -^l"- 



2kA 



- K,{x+x' —2a) 



A 



2Ka 



2na 



2na 



1 + 



A 

2na 



Jlna 



(2.6&) 
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For X, x' < 0, 



9(.,x') = ^e-"l-l + — 



Here, the denominator is 



^k{x+x') 



\' 



2Ka 



2Ka 



A 

2na 



X' 

2Ka 



2Ka 



1 



A 



1 



A' 



2K,a 



AA' 



(2.6c) 



:2.7) 



2KaJ \ 2KaJ (2fi;a)^ 

Note that in the strong couphng hmit we recover the famihar results, for example, inside 



A, A' 



oo 



g{x,x') 



sinh Kx^ sinh k(x- 



{21 



K sinh Ka 
Evidently, this Green's function vanishes at a; = and at a; = a. 

We can now calculate the force on one of the ^-function points by calculating the 
discontinuity of the stress tensor, obtained from the Green's function (j2.3|) by 

{T^'^) = Id^d"' - -g^^'d^d'A -G{x, x') 

Writing a reduced stress tensor by 

we find inside 

txx = '^{^ + 



(2.9) 



(2.10) 



4ifi:A 



-/^ 



A 



1 + 



da^t)g{x,x') 


x=x' 




f^') 


\f. ^ ^ ^ ^' ^ 2na AA' 1 

1 + 1 1 c 1 
\ 2KaJ \ 2na) (2fi;a)^ 


A' 
2K,a 




\' / \ \ 1 

X-a) 1 ^ / ^ 1 ^ \ ^2«x 

2na \ 2hiaJ 


\ 



(2.11) 



2Ka 

Let us henceforth simplify the considerations by taking the massless limit, /x = 0. Then 
the stress tensor just to the left of the point x = a is 

-1 



U 



2i 



1 + 2 



2Ka 



A 



+ 1 



2Ka 



1 1 e'^^" - 1 



(2.12a) 



From this we must subtract the stress just to the right of the point at x = a, obtained 
from (J2.6&J) . which turns out to be in the massless limit 

txx = -^, (2.126) 

x=a+ Z\ 

which just cancels the 1 in braces in ()2.12a|) . Thus the force on the point x = a due to 
the quantum fluctuations in the scalar field is given by the simple, finite expression 

_ 1 /""^ , 1 

x=a-\- 



X-'-xxl 



X-'-xxl 



Ana^ 



(iyy 



iy/X+l)iy/X' + l)ey-l 



(2.13) 



This reduces to the well-known, Liischer result jHSHnE] in the limit A, A' -^ oo. 



lim F 

A=A'^oo 



vr 



24a2 



(2.14) 
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Figure 1. Casimir force (|2.13|) between two (5- function points having strength A and 
separated by a distance a. 



and for A = A' is plotted in Fig. [T] 

Recently, Sundberg and Jaffe [HZI have used their background field method to 
calculate the Casimir force due to fermion fields between two (5-function spikes in 1 + 1 
dimension. Apart from quibbles about infinite energies, in the limit A — ;> oo they recover 
the same result as for scalar, ()2.14j) . which is as expected IHH], since in the ideal limit the 
relative factor between scalar and spinor energies is 2(1 — 2^^) in D spatial dimensions, 
i.e., 7/4 for three dimensions and 1 for one. 

We can also compute the energy density. In this simple massless case, the 
calculation appears identical, because tr^r^ = too (refiecting the conformal invariance 
of the free theory). The energy density is constant [ 1)2.111] with /i = 0] and subtracting 
from it the a-independent part that would be present if no potential were present, we 
immediate see that the total energy is E = Fa, so F = —dE /da. This result differs 
significantly from that given in Refs. yr>\\ IHHI IHH] . which is a divergent expression in 
the massless limit, not transformable into the expression found by this naive procedure. 
However, that result may be easily derived from the following expression for the total 
energy, 



E = f{dr) (r°°) = ^ f{dr){d^d'^ - V^)G{x,x') 



1 
2i 



(dr) / ^2co'gir,r) 



(2.15) 



if we integrate by parts and omit the surface term. Integrating over the Green's functions 
in the three regions, given by ()2.6aj) . ()2.66|) . and ()2.6c|) . we obtain for A = A', 

I f°° . 1 I f°° . l + 2/{y + X) 



E 



2na 



dy 



l+y/X 47ra 



dyy 



(l//A+l)V-l' 



(2.16) 
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where the first term is regarded as an irrelevant constant (A/a is constant), and the 
second is the same as that given by equation (70) of Ref. [^0] upon integration by parts. 
The origin of this discrepancy with the naive energy is the existence of a surface 
contribution to the energy. Because d^Tf^" = 0, we have, for a region V bounded by a 
surface S, 

= A /" (dr)r°° + / dS,r°\ (2.17) 



dtjv js 



Here T°* = 9°09*0, so we conclude that there is an additional contribution to the energy, 
Es = -^ f dS-VG{x,x') (2.18a) 

Zt J x'=x 

=-hr^Y.h^-y^ , ' ^2.186) 

li y_oo ^^ cla; x'=x 

where the derivative is taken at the boundaries (here x = 0, a) in the sense of the outward 
normal from the region in question. When this surface term is taken into account the 
extra terms in ()2.16p are supplied. The integrated formula ()2.15|) automatically builds in 
this surface contribution, as the implicit surface term in the integration by parts. (These 
terms are slightly unfamiliar because they do not arise in cases of Neumann or Dirichlet 
boundary conditions.) See Fulling [201 for further discussion. That the surface energy 
of an interface arises from the volume energy of a smoothed interface is demonstrated 
in Ref. [HI] , and elaborated in section 12.41 

It is interesting to consider the behavior of the force or energy for small coupling 
A. It is clear that, in fact, ()2.13|) is not analytic at A = 0. (This reflects an infrared 
divergence in the Feynman diagram calculation.) If we expand out the leading A^ term 
we are left with a divergent integral. A correct asymptotic evaluation leads to the 
behavior 

^ A^ . ., . „ A2 



^ rin2A + 7), E~-^(ln2A + 7-l), A ^ 0. (2.19) 

This behavior indeed was anticipated in earlier perturbative analyses. In Ref. [57] the 
general result was given for the Casimir energy for a D dimensional spherical 5-function 
potential (a factor of l/47r was inadvertently omitted) 

A^ F (i^) r(D - 3/2)r(l - D/2) 

Tra 2i+20[F(D/2)]2 ' ^ ' ^ 

This possesses an infrared divergence as Z^ ^ 1: 

^{D=x) ^ ^r(O), (2.21) 

47ra 

which is consistent with the nonanalytic behavior seen in (|2.19|) . 

2.2. Parallel Planes in 3 + 1 Dimensions 

It is trivial to extract the expression for the Casimir pressure between two 6 function 
planes in three spatial dimensions, where the background lies at x = and x = a. We 
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merely have to insert into the above a transverse momentum transform, 

G{x, x') = J ^e-(-*') J ^e^^(-')-,(x, x'; .), (2.22) 

where now k'^ = fi"^ + k"^ — uJ^ . Then g has exactly the same form as in ()2.fia|l - ()2.6c|l . 
The reduced stress tensor is given by, for the massless case, 

1 1 

txx = ^{dA' - K^)-g{x,x') , (2.23) 

/ 1 x=x' 

so we immediately see that the attractive pressure on the planes is given by (A = A') 

which coincides with the result given in Refs. (HDl 1^ • The leading behavior for small A 
is 

while for large A it approaches half of Casimir's result [Oj for perfectly conducting parallel 
plates. 



jTE ^ 



2 



rNJ 



A > 1. (2.256) 



480a4 
The Casimir energy per unit area again might be expected to be 

1 Pj y^ IP 



^ ~ gQn^a^ I "^y {y/X + l)2e?/ - 1 " 3 a ' ^^'^^^ 

because then P = —^£- In fact, however, it is straightforward to compute the energy 
density (T^^) is the three regions, x < 0, < x < a, and a < x, and then integrate it 
over X to obtain the energy/area, which differs from ()2.26|) because, now, there exists 
transverse momentum. We also must include the surface term ()2.18ajl . which is of 
opposite sign, and of double magnitude, to the k'^ term. The net extra term is 

2//A 



^'= .o 9 . / ^yy 



1- 



{y/\ + 1) V 



(2.27) 



487T^a^ Jo ^"^ l + y/A 
If we regard A/a as constant (so that the strength of the coupling is independent of the 
separation between the planes) we may drop the first, divergent term here as irrelevant, 
being independent of a, because y = 2Ka, and then the total energy is 

9671^ a^ Jo (i//A + l)V-l' ^ ^ 



which coincides with the massless limit of the energy first found by Bordag et al 
and given in Refs. [201 'TT] . As noted in section 12.11 this result may also readily be 
derived through use of ()2.15p . When differentiated with respect to a, ()2.28|) . with A/a 
fixed, yields the pressure (I2.24p . 

In the limit of strong coupling, we obtain 



7r2 



lim 8 = -, (2.29) 

A^oo 1440a3' ^ ' 

which is exactly one-half the energy found by Casimir for perfectly conducting plates [0]. 
Evidently, in this case, the TE modes (calculated here) and the TM modes (calculated 
in the following subsection) give equal contributions. 
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To verify this claim, we solve a similar problem with boundary conditions that the 
derivative of g is continuous at x = and a, 
d , 



g-g[x,x 



IS contmuous, 



x=Om 



but the function itself is discontinuous, 



g{x,x') 



x=a+ 



\a—g{x,x'] 



(2.30a) 



(2.306) 



and similarly at a; = 0. These boundary conditions reduce, in the limit of 
strong coupling, to Neumann boundary conditions on the planes, appropriate to 
electromagnetic TM modes: 

d 



A — i> oo : 



dx 



g{x,x') 



0. 



x=0,a 



(2.30c) 



It is completely straightforward to work out the reduced Green's function in this 
case. When both points are between the planes, < x,x' < a, 

1 J f Xna^ 



g[x,x 
+ 



1 

2k 
Xna ( _ Xna 



_^-K\x-x'\ ^ ^^i I ^w^ J 2cosh/t(x-x') 
2kA\ \ 2 ' ^ ' 



k(x+x') _|_ —K{x+x'—2a) 



2 \ 2 

while if both points are outside the planes, a < x.,x' 



(2.31a) 



g[x,x 
+ 



2k 
1 Xna 



—k\x—x'\ 



■k{x+x' —2a) 



2kA 2 
where the denominator is 



A 



1 + 



Xna 



-,2Ka 



1 - 



Xna 



Xna 



Xna 



JZna 



{2.31b) 



(2.32) 



2 / V 2 

It is easy to check that in the strong-coupling limit, the appropriate Neumann 
boundary condition ()2.30c|) is recovered. For example, in the interior region, < x,x' < 
a, 

a) 



, cosh Kx< cosh K,{x^ 

lim g[x, X ) = 

A^oo K smh Ka 



(2.33) 



Now we can compute the pressure on the plane by computing the xx component of 
the stress tensor, which is given by ()2.23p . 



txx = ^(-^^ + dxd'^)gix,x'] 



The action of derivatives on exponentials is very simple, so we find 

'ZK I Ana \ 

-K 



U 



1 
2i 



2k I Xna 

A 



2i ' 



(2.34) 

(2.35a) 
(2.356) 
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Figure 2. TE and TM Casimir pressures between 5- function planes having strength 
A and separated by a distance a. 



SO the flux of momentum deposited in tlie plane a; = a is 

in 



U 



U 



(2.36) 



.=a+ (_|_ + i) e2--l 
and then by integrating over frequency and transverse momentum we obtain the 



pressure: 



jTM 



1 



327r2a4 



dyy'^ 



1 



^ + 1 ^y 



Aj/ 



In the limit of weak coupling, this behaves as follows: 

15 



>TM 



rA^ 



(2.37) 



(2.38) 



which is to be compared with ()2.25a|l . In strong coupling, on the other hand, it has 
precisely the same limit as the TE contribution, ()2.256|) , which confirms the expectation 
given at the end of the previous subsection. Graphs of the two functions are given in 

Fig. 121 

For calibration purposes we give the Casimir pressure in practical units between 
ideal perfectly conducting parallel plates at zero temperature: 

.^_^_1.30mPa^ (2.39) 



P 



240a4 



(a/l/xm)^ 



2.4- Surface energy as hulk energy of boundary layer 

Here we show that the surface energy can be interpreted as the bulk energy of the 
boundary layer. We do this by considering a scalar field in 1 + 1 dimensions interacting 
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where 

a(x) = 



12 






2 2' 

otherwise, 



with the property that h6 = 1. The reduced Green's function satisfies 



dx'^ 



+ K, + Acr(x) 



g{x,x') = 6{x — x'). 



(2.40) 



(2.41) 



(2.42) 



This may be easily solved in the region of the slab, — 2 < x < 



2' 



g[x,x 



_]_( n'\x-x'\ ,]_ 



[K, 



/2 



K^) coshK'(x + x') 



+ {k' - nfe'^'^ cosh k'(x - x')] } . 

Here k' = ^/k^ + Xh, and 

A = 2kk' cosh k'5 + (k^ + k'^) sinh k'5. 



(2.43) 



(2.44) 



This result may also easily be derived from the multiple refiection formulas given in 
section 13. H and agrees with that given by Graham and Olum [72j • The energy of the 
slab now is obtained by integrating the energy density 



r ^ ^(.^ 



dxdx' + Xh)g 



(2.45) 



over frequency and the width of the slab. This gives the vacuum energy of the slab 
I r dK 1 



E. 



\2t 2 



2 J_^ 27r 2k'A 






S 



(2.46) 



If we now take the limit 6-^0 and /i — ^ 00 so that h6 = 1, we immediately obtain 
1 f°^ . X 



Es = — dK- 

27r Jo X + 2k 



(2.47) 



which precisely coincides with one-half the constant term in (|2.16|) . with A there replaced 
by Xa here. 

There is no surface term in the total Casimir energy as long as the slab is of finite 
width, because we may easily check that ^5'la.=^./ is continuous at the boundaries ±|. 
However, if we only consider the energy internal to the slab we encounter not only the 
energy ()2.15p but a surface term from the integration by parts. It is only this boundary 
term that gives rise to Es, ()2.47|) . in this way of proceeding. 

Further insight is provided by examining the local energy density. In this we 
follow the work of Graham and Olum [7^ 175] . However, let us proceed here with 
more generality, and consider the stress tensor with an arbitrary conformal term. 



T^- = a>a>--(7^^(aA^ 



+ A/i0")-a(9'^a'^-^'^"9- 



M^-C)2> 



(2.48) 
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in d + 2 dimensions, d being the number of transverse dimensions. Applying the 
corresponding differential operator to the Green's function ()2.43|) . introducing polar 
coordinates in the {(, k) plane, with ( = k,cos6, k = usinO, and 

(-=0 = ^. (2.49) 

we get the following form for the energy density within the slab, 

2-d-2 -{d+l)/2 roo T d , 

^" - r((/+3)/2) i ^{'^'^ - -' [(1 - ^""^ + '^^"- - ^'^ -^^'^'^ 

-(/€'-/€) V'^K^l (2.50) 



From this we can calculate the behavior of the energy density as the boundary is 
approached from the inside: 

00 Tid+l)\h l-Aaid+l)/d ,, ./. ..,,. 

~ 2'^+%('^+i)/2r((rf + 3)/2) {6-2\x\Y ' ' ' ' ^ ' 

For d = 2 for example, this agrees with the result found in Ref. [72] for a = 0: 

T'' - ^7#^^, 1^1 - -• (2.52) 

Note that, as we expect, this surface divergence vanishes for the conformal stress tensor 
[T4] . where a = d/A{d +1). (There will be subleading divergences if d > 2.) 

We can also calculate the energy density on the other side of the boundary, from 
the Green's function for x,x' < —5/2, 

. ,. 1 



g^x.x) = — 
Zk 



-k\x-x'\ k(x+x'+&)(I2 2NsinhK5 



(2.53) 



and the corresponding energy density is given by 

TOO ^ _^l_lMi±i)/^ r d«: «:'^+ii(«:'2 _ ^2y.ix+6/2) ^^^^ ^,^ ^2.54) 

2d+27r(rf+i)/2r((t/ + 3)/2) Jo /s. ' V ; 

which vanishes if the conformal value of a is used. The divergent term, as x — > —5/2, 
is just the negative of that found in ()2.5H) . This is why, when the total energy is 
computed by integrating the energy density, it is finite for rf < 2, and independent of 
a. The divergence encountered for d = 2 may be handled by renormalization of the 
interaction potential [Z2]- In the limit a.s h —>■ oo, h6 = 1, we recover the divergent 
expression ()2.47|) for d = 0, or in general 

1 /"°° A 

lim Es = ,, ^ /_,, ix;a — 7T-; r-r^ / dn k . (2.55) 

Therefore, surface divergences have an illusory character. 

For further discussion on surface divergences, see section I^?T1 
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3. Casimir Effect Between Real Materials 

3.1. The Lifshitz Formula Revisited 

As a prolegomena to the derivation of the Lifshitz formula for the Casimir force between 
parallel dielectric slabs, let us note that the results in the previous section may be easily 
derived geometrically, in terms of multiple reflections. Suppose we have translational 
invariance in the y and z directions, so in terms of reduced Green's functions, everything 
is one-dimensional. Suppose at x = and x = a we have discontinuities giving rise to 
reflection and transmission coefficients. That is, if we only had the a; = interface, the 
reduced Green's function would have the form 

g(x, x') = — (^e-"l"-"'l + re-"("+"')) , (3.1a) 

for X, x' > 0, while for x' > > x, 

^(a;,x') = — te-''(^'-^). (3.1b) 

2k 

Similarly, if we only had the interface at x = a, we would have similarly defined refiection 

and transmission coefficients r' and t'. Transmission and reflection coefficients defined 

for a wave incident from the left instead of the right will be denoted with tildes. If both 

interfaces are present, we can calculate the Green's function in the region to the right 

of the rightmost interface x,x' > a in the form 

g{x, x') = — ('e-^l^-^'l + Re-<^+^'-2^)\ , (3.2a) 

where R may be easily computed by summing multiple reflections: 

R = r' + t'e-'^Ve""'^? + t'e-^\e-''''r'e-'"'re-'"'i' + ... 

rf'f' 
= r'+ /'' . . (3.26) 

For the TE (5-function potential ()2.1|) . r = r = ~{l + 2Ka/X)^^, and t = t = 1 + r, and we 
immediately recover the result ()2.66|) . But the same formula applies to electromagnetic 
modes in a dielectric medium with two parallel interfaces, where the permittivity is 

El, X < 0, 

£:(x) = { Es, < X < a, . (3.3) 

£2, a < X. 

In that case [73] 

r = , r = , r = —r , [S.Aa) 

K,3 + Ki K2 + K3 

and 

t' = l+r', i' = l-r', (3.46) 

where k] = k"^ — 'jJ^^i- Substituting these expressions into (J3.2&J1 we obtain 

K2 - /tg 4/t2/t3 1 

K2 + /ts ^^\ - Ki ^:^^^^4^e2-3a _ l ' ^"^'^^ 
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which coincides with the formula (3.16) given in Ref. ^^. 

However, to calculate most readily the force between the slabs, we need the 
corresponding formula for the reduced Green's function between the interfaces. This 
may also be readily derived by multiple reflections: 



g[x,x) = — 



^-k\x-x'\ _|_ ~l^-K(2a-x-x') _|_ ^~l^-K{2a-x'+x) _|_ ^~/2^-K(4a-x-x') 

_|_ ^2~/2g-K{4a+x-a;') _|_ 

_^ ^g-«(^'+^) _^ ^f;'g-«(2a+a;'-x) _^ ^2~r^-K,{2a+x'+x) _^ ^2~i2^-K{4a+x' -x) _^ _ _ _ 



]_[ -n\x-x'\ 

2k\ 



+ 



f^ZiKiQi iytiyf 



2rf'coshK(x - x') + f e''(^+^') + re-'^^^+^'^^a) 



(3.6) 



Indeed, this reduces to ()2.6a|) when the appropriate reflection coefficients are inserted. 
The pressure on the planes may be computed from the discontinuity in the stress tensor, 
or 



U 



-U 



— {-K^ + d^d^>)g{x,x) 



x=x =a- 



IK 



x=x'=a+ -^e^"^" — 1 

r r' 

from which the 5-potential results ()2.12aj) and ()2.12?)|) follow immediately. For the case 
of parallel dielectric slabs the TE modes therefore contribute the following expression 
for the pressure:}:: 

""^ du f (dfc) i/«3 



P 



TE 



271 ('27r)2M±ii2«3+ae2«;3a_ X' 
The contribution from the TM modes are obtained by the replacement 



Kj ^ Kj . 



(3.^ 



(3.9) 



except in the exponentials jZHj- This gives for the force per unit area at zero temperature 



P. 



T=0 
Casimir 



-| POO pco 



(3.10) 



with the denominators here being [kj = ^yk^ + ('^ei{iC)] 

d = !!i±^!!i±^e-3« - 1, d' = 4±44±4e-3^ - 1, (3.11) 



Kg K-^ K3 K2 



1^3 - /«1 /«3 - f^2 

which correspond to the TE and TM Green's functions, respectively. This is the 
celebrated Lifshitz formula P^I12[[T^ I14J. which we shall discuss further in the following 
subsections. We merely note here that if we take the limit £1^2 — *■ 00, and set £3 = 1, we 
recover Casimir's result for the attractive force between parallel, perfectly conducting 
plates ^Ml- 

Henkel et al [70] have computed the Casimir force at short distances (~ 1 nni) from 
interactions between polaritons. Their result agrees with the Lifshitz formula with the 
plasma formula (j3.33|) employed, see Ref. [77| [TS] . 

X For the case of dielectric slabs, the propagation constant n is different on the two sides; we omit the 
term corresponding to the free propagator, however. In the energy, the omitted terms are proportional 
to the volume of each slab, and therefore correspond to the volume or bulk energy of the material. 
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3.2. The Relation to van der Waals Forces 

Now suppose the central slab consists of a tenuous medium and the surrounding medium 
is vacuum, so that the dielectric constant in the slab differs only slightly from unity, 

e-Kl. (3.12) 

Then, with a simple change of variable, 

K = Cp, (3.13) 

we can recast the Lifshitz formula ()3.10|) into the form 

P ^ --1^ rdCC'[e(C) - 1]' r %[i2p' - If + l]e-2^^^ (3.14) 

If the separation of the surfaces is large compared to the wavelength characterizing e, 
aC,c 3> 1, we can disregard the frequency dependence of the dielectric constant, and we 
find 

23(5-1)2 

For short distances, aC^^ <^ 1, the approximation is 



- / dC(e(C) - 1)^. (3.16) 

" Jo 



These formulas are identical with the well-known forces found for the complementary 
geometry in Ref. ^^. 

Now we wish to obtain these results from the sum of van der Waals forces, derivable 
from a potential of the form 

V = -^. (3.17) 

We do this by computing the energy (A/" = density of molecules) 

If we disregard the infinite self-interaction terms (analogous to dropping the volume 
energy terms in the Casimir calculation), we get jT^ ISU] 

d E _ 2nBAf^ 1 

~ ~ 5^1 ~ ~ (2 - 7)(3 - 7) a^-3 ■ ^^-^^^ 

So then, upon comparison with (|3.15|) . we set 7 = 7 and in terms of the polarizability, 

we find 

B = ^a^ (3.21) 

or, equivalently, we recover the retarded dispersion potential of Casimir and Polder |16j . 
23 a^ , , 

y = —, — 7, 3.22 

An r^ 
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whereas for short distances we recover from (j3.16p the London potential |HT] . 

V=--- dCa(C)l (3.23) 



vr r^ JO 



Recent, nonperturbative approaches to Casimir-Polder forces include that of 
Buhmann et al 



3.2.1. Force Between a Molecule and a Plate One can also calculate the force between 
a polarizable molecule, with electric polarizability ot{uj), and a dielectric slab. A simple, 
gauge-invariant way of doing this starts from the variational form |79| 1^ 

SW = - f dt5E = ~^ f{dx)S6{x)Tkkix,x), (3.24) 

where beir) = 4:TTa{uj)6{r — R), R denoting the position of the molecule. Here F is the 
electromagnetic Green's dyadic, defined by 

T{r,r') = i{E{r)E{r')). (3.25) 

In terms of the reduced Green's function, defined by fl2.22|l . then 

SE = -An / —-——a{uj)gkk{x,x-,uj,k). (3.26) 

It is easily seen how the trace of the reduced Green's function can be expressed in terms 
of the reduced TE and TM Green's functions, 



,2^TE I ^ ^TM 1 C* 1 a rpj^ 



(3.27) 



9kk=\^9 +— ^^ +-^ — 7^-7^ 
^ ee' e ox e' ox' 

For a single interface, the Green's functions to the right of a dielectric slab situated in 
the half-space x < are given by ()3.1aj) with the reflection coefficients in the vacuum 

^TE ^ ^1^^^ ^TM ^ ^-^iM ^ (328) 

K + Ki K + Ki/Si 

where k^ = k"^ + ("^ and k^ = k"^ + C^^i- I^i this way, we immediately obtain the energy 
between a dielectric slab (permittivity ei) and a polarizable molecule a distance Z from 
it: 



-1 /"OO /'OO -I 

^siab,moi = -JT^ / dC47ra(C) / dP-e- 
-Ldtt 7o Jo ^ 



-C - + i2k' + C 



(3.29) 

If the separation between the plate and the molecule is large, we expect that we may 
neglect the frequency dependence of the polarizability, a{() -^ a(0). There are then 
two simple limits. If we take ei ^ 00 we are describing a perfectly conducting plane, in 
which case we immediately obtain the result first given by Casimir and Polder J^I 

-C/metal,mol — ~T, — TTJ- [O.OU) 

On the other hand, we could consider a tenuous medium, (^i — 1) ^ 1, in which case 

_ 23 ams^-l) 

-C^dilute,mol — —-r^ ^7 ■ [6.6L 
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The latter should be, as in the previous subsection, interpretable as the sum of 
pairwise van der Waals interactions between the external molecule and the molecules 
which make up the slab, given by the Casimir- Polder interaction ()H.22j) . The net energy 
then is 

which coincides with fj3.31|) when ()3.20p is used. 

The force between a molecule and a plate has been measured by Sukenik et al. |83j . 
who actually verified the force between a molecule and two plates |S1] at the roughly 10% 
level. Recently, this result has been questioned (at about the same level of accuracy) by 
Bordag [S3] , who argued that a subtle error involving the quantization of gauge fields in 
the presence of boundaries was made by Casimir and Polder [12] and subsequent workers. 
The fact that the result can be given an unambiguous gauge-invariant derivation, and 
that it is closely related to the Lifshitz formula and the retarded dispersion van der 
Waals force suggests that this critique is invalid. (Bordag now concedes that the usual 
result is valid for "thick" plates, where the normal component of E is given by the 
surface charge density.) 

For a recent rederivation of ()3.3()j) see Hu et al jSE]- A very recent paper by Babb, 
Klimchitskaya, and Mostepanenko [87 gives a rederivation of the Casimir-Polder energy 
()3.3()|1 in the retarded limit, and finds no support for Bordag's modification. They then 
go on to discuss the dynamical polarizability and thermal corrections for real materials, 
and find substantial (35%) corrections at short distances ~ 100 nm. 

In this connection we might also mention the work of Noguez and Roman- 
Velazquez [SH] , who calculate the force between a sphere and a plate made of dissimilar 
materials in the non- retarded limit (see also van Kampen |^ and Gerlach ^^) in 
terms of multipolar interactions. They find significant deviations from the proximity 
approximation (section l3.5j) . which says that there is no difference between the force 
between a sphere made of material A and a plate made of material B and the reversed 
situation, when the separation is comparable or large compared to the radius of the 
sphere, and that under the above-mentioned A-B interchange the forces change by up 
to 6%. See also Ref. ^^. 

Ford and Sopova |93| IHl] consider Casimir forces between small metal spheres and 
dielectric (and conducting) plates, modeled by a plasma dispersion relation 

eico) = 1 - ^. (3.33) 

The electric dipole approximation used requires aup <^ 1, that is, the radius of the 
sphere a must be in the 10-100 nm range. The force is oscillatory, being alternatively 
attractive and repulsive as a function of the height Z of the atom above the plate. Thus 
levitation in the earth's gravitational field might be possible, for Z ~ 1 /im. 
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3.3. Roughness and Conductivity Corrections 

3.3.1. Roughness Corrections No real material surface is completely smooth. Even 
beyond the atomic level, there will be regions of higher and lower elevations. Insofar 
as these are plateaus large compared to the separation between the disjoint surfaces, 
the corrections can be easily incorporated by use of the proximity approximation (see 
section 13.51 below) . This is nothing other than the naively obvious statement that if 
P{a) is the force per unit area between two parallel plates separated by a distance a, 
the average force per area between rough surfaces made up of large plateaus and valleys, 
with the perpendicular distance between two adjacent points on the two surfaces in terms 
of transverse coordinates {x,y) being a{x,y), is 

P=jJdxdyP{a{x,y)). (3.34) 

In Ref. jnS], for example, an equivalent expression is used directly with data obtained by 
topography of the surfaces using an atomic force microscope. Traditionally, a stochastic 
estimate has been used. Let the separations a be distributed around the mean ao 
according to a Gaussian, with the probability of finding separation a being given by 

p(a) = — i-e-^"-"")'/^^'^)'. (3.35) 

^/7^oa 

We will assume 5a <^ oq. Then, (a) = ao, ((a — ao)^) = ^{Sa^, and in general 

/•OO -I /"CO 

{a")= daaXa) = ^.— / dae-'^'/^^'^Va + ao)° 
Jo V^oa J_^ 



% 



a{a-l)l{5af a{a - l){a - 2){a - 3) 3 {da^ 

2 2^r^ 4! 4^r^' 



(3.36) 

The force between a sphere and a plate depends on the closest distance d between them 
like d~^, see fl3.78|) below, so the stochastic estimate for the roughness correction in that 
case, in terms of the mean-square fluctuation amplitude A = Sa/y/2, is 



sph— pi, rough -f^sph— pi 



A^ 


2 


(A^ 


4 




+ 45 




+ ... 


d, 




Vdy 





1 + 6- +45- +... . (3.37) 



A much more detailed discussion may be found in Ref. jTUI. It must be appreciated that 
the approximate treatment based on the proximity approximation is invalid for short 
wavelength deformations 



3.3.2. Finite Conductivity Another interesting result, important for the recent 
experiments ^T7\ IHHl IHni HOOj - is the correction for an imperfect conductor, where for 
frequencies above the infrared, an adequate representation for the dielectric constant is 
[75] that given by the plasma model ()3.33|1 where the plasma frequency is, in Gaussian 
units 



Ane'^N 



^y = , (3.38) 
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where e and m are the charge and mass of the electron, and A^ is the number density 
of free electrons in the conductor. A simple calculation shows, at zero temperature 

mniizn], 



p 



TT 



240a4 



1 



m\ 1/2 

N 



3y^ ea 
If we define a penetration parameter, or skin depth, hy 6 = l/wp, 
per area for parallel plates out to fourth order as |1U21 EOJ I1U3| ITU] 



(3.39) 

we can write the force 



P 



TC 



240a4 



■ 16 6 ^5'^ 640 
1 - h24— - 



vr 



\ 5^ 2800 



1 



leSvr^ 



i4 



(3.40) 



3 a a? 7 V 210y a^ 9 V 7350 J a-* 

while using the proximity force theorem (see section l375|l . to convert pressures between 
parallel plates to forces between a lens of radius R and a plate, 

27ri? 



^, 



n-l 



n — 1 



-aP„ 



for a term in the pressure going like Pn (x. a 
and a plate separated by a distance d is 



(3.41) 
the force between a spherical surface 



T 



TC^R 



5 72 5^ 
d^ 5 d^ 



320 



1 



71 

210 



t/3 



400 



1 



1637r^ 



rf4 



(3.42) 



360t/3 [" "t/ ■ 5 (^2 7 \^- 210y t/3 ■ 3 \^- 7350 

Lambrecht, Jaekel, and Reynaud jl04j analyzed the Casimir force between mirrors 
with arbitrary frequency-dependent reflectivity, and found that it is always smaller than 
that between perfect reflectors. 

We might also mention here the interesting suggestion that repulsive Casimir forces 
might exist |105j between parallel plates. This harks back to an old suggestion of Boyer 
|l()(ij , that repulsion will occur between two plates, one of which is a perfect electrical 
conductor, e — > cxd, and the other a perfect magnetic conductor, fi —^ 00, 

P=-~ (3.43) 

8 240 a^ ^ ' 

However, it appears that it will prove very difficult to observe such effects in the 
laboratory |I07j . Klich |l()8j now seems to agree with this assessment. 



3.4- Thermal Corrections 

The discussion in this subsection is adapted from that in Refs. |1()9| IllOj . We begin 
by reviewing how temperature effects are incorporated into the expression for the 
force between parallel dielectric (or conducting) plates separated by a distance a. To 
obtain the finite temperature Casimir force from the zero-temperature expression, one 
conventionally makes the following substitution in the imaginary frequency, 

C - C^ = ^, (3.44a) 

and replaces the integral over frequencies by a sum, 

dC 1 ^ 

> 27r ^ /? 



E- 



(3.446) 
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This reflects the requirement that thermal Green's functions be periodic in imaginary 
time with period f3 |lllj . Suppose we write the finite-temperature pressure as [for the 
explicit form, see ()3.1()|1 and ()3.58j) below] 

oo 

P^ = J2'fm, (3.45) 

m=0 

where the prime on the summation sign means that the m = term is counted with 
half weight. To get the low temperature limit, one can use the Euler-Maclaurin (EM) 
sum formula, 



oo 



/•oo -| R 

J2 fik) = / m dk + -m - J2 77^/''-'\o), (3.46) 

fc=0 "^0 ^ 9=1 ^^^>- 

where -B„ is the nth Bernoulli number. This means here, with half-weight for the m = 
term, 

/oo °° TD 

fim)dm-J2^f^^'-'\0). (3.47) 

fc=i ^ ^' 

It is noteworthy that the terms involving /(O) cancel in ()3.47|) . The reason for this is 
that the EM formula equates an integral to its trapezoidal-rule approximation plus a 
series of corrections; thus the 1/2 for m = in ()3.45|1 is built in automatically. For 
perfectly conducting plates separated by vacuum [see the A ^ cxo limit of ()2.24|) or 
flTTfjl . or the £1,2 -^ oo limit of (jTTTl|l with £3 = 1] 

f(x) = -^ K^dK—^ . (3.48) 

Of course, the integral in ()3.47j) is just the inverse of the finite-temperature prescription 
()3.446|) . and gives the zero-temperature result. The only nonzero odd derivative 
occurring is 

no) = -^, (3.49) 

which gives a Stefan's law type of term, seen in (j3.53|) below. 

The problem is that the EM formula only applies if f{fn) is continuous. If we follow 
the argument of Ref. |!?K t l!?H t li ^lll2j . and take the ei,2 -^ 00 limit of ()3.1()j) at the end§ 
(ei,2 are the permittivities of the two parallel dielectric slabs), this is not the case, and 
for the TE mode 

/o = 0, (3.50a) 

C(3) 2'Kam 

Then we have to modify the argument as follows: 

00 00 00 ^ 

^^ = E '/- = E /- = E 'fm - 2/0, (3-51) 

m=0 m=l m=0 

§ This is contrary to the "Schwinger" prescription advocated in Refs. IZniEHl, in which the perfect- 
conductor Umit is taken before the zero-mode is extracted. 
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where /„ is defined by continuity, 

fm={ ,. ''^\. "" ^ °' (3.52) 

— , m = 0. 




Then by using the EM formula, 







240a4 



' + TV/3 



C(3) n^ 1 
l7r/5a3 45 /^^ 

C(3) 



16 / a^ 



ina^ 



aT < 1. (3.53) 



The same result for the low-temperature limit is extracted through use of the Poisson 
sum formula, as, for example, discussed in Ref. ^^. Let us refer to these results, with 
the TE zero mode excluded, as the modified ideal metal model (MIM). The conventional 
result for an ideal metal (IM), obtained first by Lifshitz fU E] and by Sauer [21] and 
Mehra [22]; is given by (j3.53p with the linear term in T omitted. 

Exclusion of the TE zero mode will reduce the linear dependence at high 
temperature by a factor of two, 

^"--S^- ^"M~-|§r, aT»l, (3.54) 

but this is not observable by present experiments. The observable consequence, however, 
is that it adds a linear term at low temperature, which is given in ()3.53|1 . up to 
exponentially small corrections |^ . 

There are apparently two serious problems with the result ()3.53|1 : 

• It would seem to be ruled out by experiment. The ratio of the linear term to the 
T = term is 

A = ^^^^aT = 1.16aT, (3.55a) 

or putting in the numbers (300 K = (38.7)~^ eV, he = 197 MeV fm) 

A = 0.15 (—^] (^-] , (3.556) 

V300K J Vl/imy ' ^ ' 

or as Klimchitskaya observed |113j . there is a 15% effect at room temperature at 
a separation of one micron. One would have expected this to have been been seen 
by Lamoreaux |114j : his experiment was reported to be in agreement with the 
conventional theoretical prediction at the level of 5%. (Lamoreaux |115j is now 
proposing a new experiment to resolve this issue.) 

• Another serious problem is the apparent thermodynamic inconsistency. A linear 
term in the force implies a linear term in the free energy (per unit area), 

^ = ^0 + 4^^' «^ ^ 1' (3-56) 

lovra^ 

which implies a nonzero contribution to the entropy/ area at zero temperature: 
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Taken at face value, this statement appears to be incorrect. We will discuss this problem 
more closely in section l3.4.3l and will find that although a linear temperature dependence 
will occur in the free energy at room temperature, the entropy will go to zero as the 
temperature goes to zero. The point is that the free energy F for a finite e always will 
have a zero slope at T = 0, thus ensuring that 5 = at T = 0. The apparent confiict 
with ()3.57p or ()3.53|) is due to the fact that the curvature of F{T) near T = becomes 
infinite when e —>■ oo. So ()3.56|) and ()3.57|) . corresponding to the modified ideal metal 
model, describe real metals approximately only for low, but not zero temperature - See 
the following. 

3.4-1 ■ Lifshitz formula at nonzero temperature The Casimir surface pressure at finite 
temperature P^ between two dielectric plates separated by a distance a can be obtained 
from the Lifshitz formula ()3.10|) by the prescription ()3.445|) || 



oo 



P^ = -^J2' ^"^'^ [(^™ ^''^^ - 1)"' + {B^'^''"' - 1)"'] • (3-58) 

The relation between k and the transverse wave vector k± is k^ = k'j_ + (^, where 
(rn = ^nm/iS. Furthermore, the squared refiection coefficients are 



Am — 

\ep+sj 


Bm = 


f s — p 
\s+p 


S^ =£-!+/, 


K 
P=. , 





2 



(3.59a) 
(3.596) 

with e{i(rn) being the permittivity. Here, the first term in the square brackets in (|3.58|) 
corresponds to TM modes, the second to TE modes. Note that whenever e is constant. 
Am and Bm depend on m and k, only in the combination p, 

Amin) = A{p), Bm{K) = B{p). (3.60) 

The free energy F per unit area can be obtained from ()3.58p by integration with 
respect to a since P^ = —dF/da. We get |117j 

(3F = —Y' Kdfi:[ln(l-A™)+ln(l-A™)], (3.61a) 



=0 -^C, 



where 



A™ = A^e-2'^^ A^^ = Bme-^""". (3.616) 

From thermodynamics the entropy S and internal energy U (both per unit area) 
are related io F hy F = U — TS, implying 

S=-^, and thus U = ^^. (3.62) 

II A rederivation of the Casimir force between dissipative metallic mirrors at nonzero temperature has 
been given by Reynaud, Lambrecht, and Genet |47j . They obtain formulas, generalizing those at zero 
temperature |116| . for the force valid even if the smoothness condition necessary for the derivation of 
the Lifshitz formula is not satisfied due to the failure of the Poisson summation formula. 
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As mentioned above the behaviour of S" as T — i> has been disputed, especially for 
metals where e — > oo. We now see the mathematical root of the problem: The quantities 
Am = Bm ^ 1 in the e — *> oo limit except that Bq = Q for any finite e. So the question 
has been whether Sq = or i?o = 1 or something in between should be used in this 
limit as results will differ for finite T, producing, as we saw above, a difference in the 
force linear in T. The corresponding difference in entropy will thus be nonzero. Such 
a difference would lead to a violation of the third law of thermodynamics, which states 
that the entropy of a system with a nondegenerate ground state should be zero at T = 0. 
Inclusion of the interaction between the plates at different separations cannot change 
this general property. We will show that this discrepancy vanishes when the limit £ ^ cxd 
is considered carefully. 

3.J1..2. Gold as a numerical example Let us go back to ()3.58p for the surface pressure, 
making use of the best available experimental results for e{iC,) as input when calculating 
the coefficients A^, and B^,- We choose gold as an example. Useful information about 
the real and imaginary parts, n' and n", of the complex permittivity n = n' + in", versus 
the real frequency uj, is given in Palik's book |118j and similar sources. The range of 
photon energies given in Ref. |118j is from 0.1 eV to 10^ eV. (The conversion factor 

1 eV = 1.519 X 10^^ rad/s (3.63) 

is useful to have in mind.) When n' and n" are known the permittivity e{i() along the 
positive imaginary frequency axis, which is a real quantity, can be calculated by means 
of the Kramers-Kronig relations. 

Figure El shows how £:(iC) varies with ( over seven decades, C G [10"'^^, 10^^] rad/s. 
The curve was given in Refs. fTll I119J . and is reproduced here for convenience. (We 
are grateful to A. Lambrecht and S. Reynaud for having given us the results of their 
accurate calculations.) At low photon energies, below about 1 eV, the data are well 
described by the Drude model, 

where ujp is the plasma frequency ()3.38p and u the relaxation frequency. (Usually, z/ is 
taken to be a constant, equal to its room-temperature value, but see below.) The values 
appropriate for gold at room temperature are [771 I119J 

ujp = 9.0eV, z/ = 35 meV. (3.65) 

The curve in Fig. |21 shows a monotonic decrease of ^(iC) with increasing (, as 
any permittivity as a function of imaginary frequency has to follow according to 
thermodynamical requirements. The two dashed curves in the figure show, for 
comparison, how e{i(,T) varies with frequency if we accept the Drude model for all 
frequencies, and include the temperature dependence of the relaxation frequency with 
T as a parameter. (The latter is given in Fig. ^ according to the Bloch-Griineisen 
formula |120j , which, however, does not take into account the physical fact that because 
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^ (rad/s) 

Figure 3. Solid line: Permittivity e{i() as function of imaginary frequency ( for 
gold. The curve is calculated on the basis of experimental data. Courtesy of Astrid 
Lambrecht and Serge Reynaud. Dashed lines: e(iC) versus ^ with T as parameter, 
based upon the temperature dependent Drude model; cf. Appendix D of Ref. |l()9j . 
The upper curve is for T = 10 K; the lower is for T = 300 K, which for energies below 
1 eV (1.5 X 10^^ rad/s) nicely fits the experimental data. Both curves are below the 
experimental one for ^ > 2 x 10^^ rad/s. 



of impurities, no actual conductor has zero resistivity at zero temperature |121j . See 
Appendix D of Ref. |l()9j .) For T = 300 K, the Drude curve is seen to be good for all 
frequencies up to C ~ 2 x 10^^ rad/s; for higher ( it gives too low values of e. Both 
Drude curves, for T = 10 K and T = 300 K, are seen to give the same values when 
C > 3 X 10^^ rad/s. 

As experiments are usually made at room temperature for various gap widths, we 
show in Fig. El how the surface force density for gold varies with a, at T = 300 K. The 
linear slope seen for a > 4/im is nearly that predicted by ()3.54j) for high temperatures 
when the TE zero-mode is excluded (modified ideal metal), which gives a slope of 
2.0 X 10-2^ NmVAim. (This is in spite of the fact that aT = 0.5 at a = 4 /xm.) The 
linear region between 1 and 2 /xm corresponds roughly to that in ()3.53|1 (intermediate 
temperatures). Also shown is the prediction of the temperature dependent Drude model, 
when T = 300 K. The differences are seen to be very small. Since the Drude values for 
the permittivity are lower than the empirical ones at high frequencies, as seen in Fig. 01 
we expect the predicted Drude forces to be slightly weaker than those based upon the 
empirical permittivities. This expectation is borne out in Fig. El the differences being 
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300 



Figure 4. Temperature dependence of the relaxation frequency for gold based on the 
Bloch-Griineisen formula |12()| . 



large enough to be slightly visible at short distances, as we would expect since the plasma 
nature of the material becomes more pronounced for small distances. Note that the 
temperature dependence of the permittivity is irrelevant here because the temperature 
is fixed. 

It is of interest to check the magnitude of the dispersive effect in these cases. 
We have therefore made a separate calculation of the expression ()3.58j) when e is 



taken to be constant. Figure IHl shows how the force varies with aT in cases when 
e G {100, 1000, 10000, oo} are inserted in the expressions for Am and B^ in ()3.59aj) . 

It is seen from the figure that the first three curves asymptotically approach the 
£ = oo curve, when e increases, as we would expect. Again, we emphasize that the 
dispersive curve for gold is calculated using the available room-temperature data for e{i() 
from Fig. El In the nondispersive case, there is of course no permittivity temperature 
problem since e is taken to be the same for all T. 

There are several points worth noticing from Fig. IHl 

(i) The curves have a horizontal slope at T = 0. For finite e this property is clearly 
visible on the curves. This has to be so on physical grounds: If the pressure had 
a linear dependence on T for small T so would the free energy F, in contradiction 
with the requirement that the entropy S = —dF/dT has to go to zero as T ^ 0. 
For the gold data the initial horizontal slope is not resolvable on the scale of this 
graph, but see the discussion in section 15.4.31 
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T = 300 K 




Figure 5. Surface pressure for gold, multiplied with a^, versus a when T — 300 K. 
Input data for e(iC) are taken from Fig. |21 



(ii) The curves show that the magnitude of the force diminishes with increasing 
T (for a fixed a), in a certain temperature interval up to aT ^ 0.3. This 
perhaps counterintuitive effect is thus clear from the nondispersive curves. This is 
qualitatively similar to the behavior seen in Fig. IHlfor fixed T, where the minimum 
occurs for aT ~ 0.4. 

(iii) It is seen that the curve for e = const. = 1000 gives a reasonably good 
approximation to the real dispersive curve for gold when a = 1 /im; the deviations 
are less than about 5% except for the lowest values of aT {aT < 0.1). This 
fact makes our neglect of the temperature dependence of e{iQ appear physically 
reasonable; the various curves turn out to be rather insensitive with respect to 
variations in the input values of ^(^C)- 

(iv) Also, it can be remarked that i?o = is required when e is finite. Otherwise the 
curves in Fig. IHl and thus the free energy, would have a finite slope at T = which 
again would imply a finite entropy contribution at T = in violation with the third 
law of thermodynamics. 



3.4-3. Behavior of the Free Energy at Low Temperature The low temperature 
correction is dominated by low frequencies,^ where the Drude formula is extremely 

% This statement is in the context of using of the Euler-Maclaurin summation formula to evaluate 
(|3.58|) . for example. 
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Figure 6. Nondispersive theory: Surface pressure for e E {100, 1000, 10000, oo}. 
For fow values of aT the latter coincides with the expression H3.53|l . Also shown for 
comparison is the dispersive result for gold, where experimental input data for e{\Q 
are taken from Fig. |2| Gap width is a = 1 /im. The constraint a = 1 /im applies only 
to the dispersive case, since otherwise a'^P^ is a function of aT only. Note that room 
temperature (300 K) corresponds to aT = 0.13. 



accurate. Using this fact, we have performed analytic and numerical calculations which 
show that the free energy has a quadratic low-temperature dependence, independent of 
the plate separation: 

2 



4oZ/ 



l) = Fo + T2(19eV), T< 



CUpfl^ 



20 mK 



(3.66) 



where we have put in the numbers for gold, ()3.65|) . (the temperature restriction refers 
to a 1 /im plate separation) rather than the naive extrapolation (j3.56|) 



F = Fo + T 



C(3) 



Fo + 



T 



-0.30. 



(3.67) 



We see from Fig. [7| that the value in (j3.67p indeed results if one extrapolates the 
approximately linear curve there for (a > 0.25 to zero, following the argument given in 
(|3.5ip . However, we see that the free energy smoothly changes to the quadratic behavior 
exhibited in ()3.66|) . Of course, the turn-over will be much sharper if we replace the room- 
temperature relaxation frequency z/(300 K) by the positive value at zero temperature, 
due to elastic scattering from defects or impurities. 

Results consistent with these have been reported by Sernelius and Bostrom |122j . 
In particular they show that one cannot ignore the constant value z/(0 K), so there is 
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Figure 7. The behavior of the free energy for low frequencies, in the Drude model, 
with parameters suitable for gold, and a plate separation of a = 1 fiui. Here 
^^^ ~ 2t^ Sm=o'/(Cm)- Here, we have used the room temperature value of the 
relaxation parameter. 



no relevant temperature dependence of the relaxation parameter. Although there is a 
region of negative entropy, the Nernst heat theorem is not violated, but rather S -^ 
as T ^ if one goes to sufficiently low temperature, in contradiction to Refs. [1221 ESI • 



3.4-4- Surface impedance form of reflection coefficient It has been proposed that 
the resolution to the temperature problem for the Casimir effect is that the surface 
impedance form of the reflection coefficients should be used in the Lifshitz formula 
|1241 11251 11261 IT77J ■ rather than that based on the bulk permittivity. Here we show that 
the two approaches are in fact equivalent, and that the former must include transverse 
momentum dependence. 

For the TE modes, the reflection coefficient is given by ()H.4aj) [ 7H] 



„TE 



where 



Kn 



kiz + hz ' 



uS^e — k\ 



iVCRiCFT] 



%Kr, 



(3.68) 



(3.69) 



with K^ = Kg = A;^ + C^) and the subscripts 1 and 2 refer to the metal and the vacuum 
regions, respectively. Now from Maxwell's equations outside sources we easily derive 
just inside the metal (the tangential components, designated by _L, of E and B are 
continuous across the interface) 



i^iz^j. ■ -B_L — Vjje I 1 






fci 



(n X Ei_) = 0, 



(3.70 a) 
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-iki^k^ ■ (n X Ei_) - iuk^ -3^ = 0. (3.706) 
Here n is the normal to the interface. Now the surface impedance is defined by 

E± = Z{uj, k^)Bx_ X n, n x E± = Z{uj, k^)B^. (3.71) 
So ehminating B± using this definition we find two equations: 

ki. = -|, (3.72) 

kl^ = cu^e _ kl, (3.73) 

the latter being the expected dispersion relation ()3.69|) . Substituting this into the 
expression for the reflection coefficient ()3.(i8j) we find 



TE C + Zk 1 + Zp K 

Q — Zk 1 — Z-p Q 

which apart from (relative) signs (presumably just a different convention choice) 
coincides with that given in Geyer et al |124j or Bezerra et al |128j . See also 
Refs. |129| I13()j . The first discussion of the Lifshitz formula in this approach was given 
in Ref. [TO^ . 

However, it is crucial to note that the "surface impedance" so defined depends on 
the transverse momentum, 

Z = ^ ^, (3.75) 

and so r"^^ — ;> as C ~^ just as in the dielectric constant formulation. Of course, we 
have exactly the same result for the energy as before, since this is nothing but a slight 
change of notation, as noted in Ref. |13H IHH] . 

It is therefore incorrect to assume that Z is only a function of frequency, not of 
transverse momentum, and to use the normal and anomalous skin effect formulas derived 
for real waves impinging on imperfect conductors."*" In the above-cited references, this 
necessary dependence was not included. (For further comments on the insufficiency of 
the argument in Ref. |124j see Ref. |132j . 

How does the usual argument go? The normal component of the wavevector in a 
conductor is given by 



k.. 






1/2 



Vi4:TTiua, cij ^ 0, (3.76) 



from which the usual normal skin effect formula follows immediately. 



Z{u) = -{l-i)J^. (3.77) 

However, the limit in ()3.76p here consists in omitting two "small" terms: u'^e (which is 
legitimate) and k'j_ < u"^. Here this last is not valid because in going to finite temperature 

"*" Of course, in general, the permittivity will be a function both of the frequency and the transverse 
momentum, e{uj, k±), but we believe the latter dependence is not significant for separations larger than 
hc/up = 0.02 /im. 
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we have severed the connection between u ^ i( and k±; the latter is in no sense ignorable 
as we take C ^ to determine the low temperature dependence. This is the same error 
to which we refer in Ref. J1U.9J • (This k± dependence still seems to be ignored in a 
recent reanalysis by Torgerson and Lamoreaux |133j (see also Ref. |134j ) who argue 
that low frequencies of order of the inverse transverse size of the plates dominate the 
low temperature behavior so that a linear term in the temperature does not appear. 
This seems unlikely since the zero-temperature dependence is extracted by an analytic 
continuation procedure.) 

Not only do Mostepanenko, Klimchitskaya, et al |129| I124J ignore transverse 
momentum dependence, but they apparently do not use the correct values of the 
frequency in their evaluation of the surface impedance. They use the impedance 
appropriate to the domain of infrared optics, thereby extrapolating the surface 
impedance at what they consider a characteristic frequency ~ l/2a rather than using the 
actual zero frequency value |126j . This seems to be a completely ad hoc prescription, 
as opposed to the procedure advocated in Brevik et al |l()9j . which uses the actual 
electrical properties of the materials. 

A beginning of a general discussion of nonlocal effects, including the anomalous 
skin effect, in Casimir phenomena has recently been given by Esquivel and Svetovoy 
|135j . There they argue that the Leontovich approach |136t I137J advocated by |129| I124J 
only applies to normal incidence, which is why the surface impedances only depend on 
frequency. In fact, this is incorrect in general, and if only local functions are used for the 
permittivity, that is £ = ^(a;), the dependence for the TE surface impedance given above 
is reproduced. For propagating waves the Leontovich approximation is appropriate, but 
not for the evanescent fields relevant to the Casimir effect, where k±/uj > 1 occur. They 
do not calculate temperature effects; the nonlocal anomalous skin effect for uj < oOp that 
they compute gives a correction to the Casimir force of order 0.5%, but other nonlocal 
effects, such as plasmon excitations, could be more significant |138[ll39j . 

3.5. Beyond the Proximity Approximation 

As we will discuss in the next section, to avoid problems of parallelism, most recent 
experiments to measure the force between conductors have not been made between 
parallel plates, but between a plate and a spherical surface, or between crossed cylinders. 
The Lifshitz and Casimir formulas do not apply to these situations. However, in the 
1930s, it was recognized that if the separation between the sphere and the plate is very 
small compared to the radius of curvature of the sphere, the latter force may be derived 
from the force for the parallel plate configuration. This result is usually called the 
Proximity Force Theorem |140j . which here says that the attractive force JF between a 
sphere of radius R and a flat surface is simply the circumference of the sphere times the 
energy per unit area for parallel plates, or, from ()2.29|) . 

^=2.fl£M = -|lf|, U»,, (3.78) 
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where d is the distance between the plate and the sphere at the point of closest approach, 
and R is the radius of curvature of the sphere at that point. (The exact shape of the 
"sphere" is not relevant in the strict approximation R^ d.) The proof of ()3.78j) is quite 
simple, li R^ d, each element of the sphere may be regarded as parallel to the plane, 
so the potential energy of the sphere is 

V{d)= 2nRsm9Rd9E{d + R{l-cose)) = 2nR dx£{d + R-x). (3.79) 

Jo J-R 

To obtain the force between the sphere and the plate, we differentiate with respect to 
d: 

dv r^ d 

T= -—- = 2txR dx—S(d + R-x) 
ad J_ji ox 

= 2nR[£{d)-£{d + 2R)]^2nR£{d), d <^ R, (3.80) 

provided that £{a) falls off with a. This result was already given in Refs. |14Hll42(ll43j . 
The proximity theorem itself dates back to a paper by Derjaguin in 1934 |144[I145] . 

Let us apply this theorem to the MIM model ()3.53|) for the force between parallel 
plates at low temperature. The corresponding free energy is 

- - ""' ' ^aT' - ^T' + ^^T, (3.81) 



720a3 45 27r levra^' 

where the term constant in a is determined by the high-temperature limit ()3.54|) - see 
Ref. ^ni; P- 56. This free energy is to be used in the proximity force theorem, with the 
result for the force between a sphere and a plate [HHl 11141 I146t [TUl 






(3.82) 



360 rf3 

The terms linear in T would not be present in the IM model. At room temperature, 300 
K, and at 1 /im separation, the successive terms correspond to corrections of —0.46%, 
+3.1%, and —23%, respectively. This model, of course, does not begin to reflect the true 
temperature dependence, discussed for parallel plates above. A full discussion of the 
temperature dependence for the force between a spherical lens and a plate will appear 
elsewhere. 

Emig has recently presented exact results for Casimir forces between periodically 
deformed surfaces |147l I148| 1149] . In the latest paper, the authors calculate the force 
between a fiat plate and one with a rectangular (square) corrugation, of amplitude Aa. 
This was probed experimentally by Roy and Mohideen |150j , with clear deviations from 
the proximity approximation. (See also Refs. |15H I152J for measurements of the so- 
called "lateral Casimir effect.") For short wavelength corrugations for either TE or TM 
modes one gets 

7r2 1 7r2 / 4Aa\ 

P = 7 ^^ ~ J 1 + , 3.83 

480 (a - Aa)4 4800^ \ a J ' ^ ' 

while for long wavelength corrugations 

480 2 V(a-Aa)4^ (a + Aa)V ' ^-^-^ ^ 
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which is as expected from the proximity approximation. For intermediate wavelength 
corrugations numerical results are given. The force approaches that given by the 
proximity approximation for large A like Aa/A, as compared to (Aa/A)^ for sinusoidal 
corrugations, due to the sharp edges. These behaviors can be understood from the ray 
optics approach of Jaffe and Scardicchio |2S1 discussed in the following subsection. The 
relative contributions of the TE and TM modes vary with the wavelength and the shape 
of the corrugation, the ratio of the modes approaching unity as a/Aa tends to 1 or oo. 
Insofar as first approximations to these interactions were extracted through use of the 
proximity force theorem, these results shed valuable light on how to move beyond that 
approximation. 

3.5.1. Optical Paths A very interesting strategy for moving beyond the proximity 
approximation has been suggested by Jaffe and Scardicchio [2H]. This is related to the 
semiclassical closed orbit approach advocated by Schaden and Spruch |153t I154| 1155] 
and earlier by Gutzweiler |156| I157J , and also to that of Balian and Bloch jl58t I159| I160J . 
Fulling has also recently proposed similar ideas |16Hll(i2J . 

In the simplest context, that of parallel plates, the approach is, of course, exact, 
and is precisely what we wrote down in ()3.fi|) . We simply compute the energy using 
(jTTH|l with 

gir,r) = l^^g{x,x), (3.85) 

where g{x,x') is given by ()3.fij) . Rather than carry out the sum as given there, let us 
sum the terms with even and odd numbers of reflections separately. The former give, 
when the zero reflection term is omitted, 

ge.en{x,x) = ^2 [r'TB'^'^^ + {r'rfe-''''' + ...]= ;^(coth«:a - 1), (3.86) 

where in the last step we have inserted the values for the reflection amplitudes 
appropriate to Dirichlet boundaries, r = f' = — 1. When this is inserted into the 
expression for the energy we obtain rather immediately the usual result for the Casimir 
energy between Dirichlet plates: 

S = -^ / du—^- = -—. (3.87) 

967r2a3 Jq e" - 1 14400^ ^ ^ 

Keeping only the first term in the sum (2 reflections) gives 

^^'^ = -77^' (3.88) 

which is in magnitude only 7.6% low, while keeping 2 plus 4 reflections give an error of 

1.8%: 

The odd reflections give a term in g{x,x) which depends on x: 

9om{x, X) = -i- {e-'^^ + e^'^(-")) ^^^. (3.90) 
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when this is integrated over x, the a dependence of this term disappears, so this gives 
rise to an irrelevant constant in the energy. Keeping it and the zero-reflection term gives 
the expression for the total energy as obtained directly from (jH.fij) 

1' 



So + Seven + ^odd = ~T^~ / ^V V^ [ cothy ~ -] ■ (3.91) 



Jaffe and Scardicchio j23] use this method to estimate the force between a sphere 
an a plate. The results disagree with the proximity approximation when d/R is bigger 
than a few percent, but agrees with an exact numerical calculation |163j . described in 
the following subsection, up to d/R ~ 0.1, where the proximity theorem fails badly. 

3.5.2. Worldline Approach to the Casimir Energy Gies, Moyaerts, and Langfeld 
[163, I164J have developed a numerical technique for extracting Casimir energies in 
nontrivial geometries, such as between a sphere and a plate. It is based on the string- 
inspired worldline approach. They consider, like Graham et al |165t IHTl I3U] a scalar 
field in a smooth background potential like ()2.40|) . The worldline representation of 
the effective action is obtained by introducing a proper time representation of the 
functional logarithm with ultraviolet regularization, doing the trace in configuration 
space, and interpreting the matrix element there as a Feynman path integral over all 
worldlines x(r). Field theoretic divergences can thus be handled. Other divergences 
arise from the potential itself, when it approaches some idealized limit, which may not 
be removed in a physically meaningful way and may or may not contribute to physical 
observables. The expectation value is evaluated by the "loop-cloud" method, using 
techniques from statistical mechanics. Although in the "sharp" and "strong" limits in 
the sense of Graham et al |165t 1^ IHU] divergences occur in the theory, a finite force 
between rigid bodies can be obtained. The general result for (5-function planes, discussed 
in section 1221 is reproduced numerically, and then the sphere-plate system is considered. 
The numerical results, for d/R from 10~^ to 10, agree closely with the geometric mean* 
of the plate-based and the sphere-based proximity force approximation (deviation from 
either becomes sizable for d/R > 0.02). Note that electromagnetic fluctuations (e.g., 
TM modes) have not been considered in this approach. 

3. 6. Status of the Experimental Measurements on Casimir Forces 

Attempts to measure the Casimir effect between solid bodies date back to the 
middle 1950s. The early measurements were, not surprisingly, somewhat inconclusive 
[1121 ESI CHHl CSni IHni imi III2I III3 imi nisi- The Llfshltz theory (nmUl . for zero 
temperature, was, however, confirmed accurately in the experiment of Sabisky and 

* The geometric mean version of the proximity force approximation, which coincides with the 
semiclassical periodic orbit method of Schaden and Spruch |153[ 11541 I155| . has been found to be 
the most accurate also for concentric cyhndrical shells, the Casimir energy for which was calculated by 
Mazzitelh et al [TCS1ITB7I. 
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Anderson in 1973 [IH]. So there could be no serious doubt of the reahty of zero-point 
fluctuation forces. For a review of the earlier experiments, see Refs. P^ll76j . 

New technological developments allowed for dramatic improvements in experimen- 
tal techniques in recent years, and thereby permitted nearly direct confirmation of the 
Casimir force between parallel conductors. First, in 1997 Lamoreaux used a electrome- 
chanical system based on a torsion pendulum to measure the force between a conducting 
plate and a sphere |114[ll4b'] . as given by the proximity force theorem ()3.78p . Lamore- 
aux |114[I146] claimed agreement with this theoretical value at the 5% level, although it 
seems that finite conductivity was not included correctly, nor were roughness corrections 
incorporated |177j . Further, Lambrecht and Reynaud [ZZj analyzed the effect of conduc- 
tivity and found discrepancies with Lamoreaux |178j . and therefore stated that it was 
too early to claim agreement between theory and experiment. See also Refs. |119|I179] . 

An improved experimental measurement was reported in 1998 by Mohideen 
and Roy jHZI, based on the use of an atomic force microscope. They included 
finite conductivity, roughness, and conventional temperature corrections, although no 
evidence for latter has been claimed. Spectacular agreement with theory at the 1% 
level was attained. Improvements were subsequently reported [031001 • (The nontrivial 
effects of corrugations in the surface were examined in Ref. [ISUI I15H I152J .) Erdeth 
|18Uj measured the Casimir forces between crossed cylinders at separations of 20-100 
nm. The highest precision was achieved with very smooth, gold-plated surfaces. Rather 
complete analyses of the roughness, conductivity, and temperature corrections to the 
Lamoreaux and Mohideen experiments have been published J181L I182[ EO] • 

More recently, a new measurement of the Casimir force ()3.78p was presented by 
a group at Bell Labs |183[ I184J . using a micromachined torsional device, a micro- 
electromechanical system or MEMS, by which they measured the attraction between 
a polysilicon plate and a spherical metallic surface. Both surfaces were plated 
with a 200 nm film of gold. The authors included finite conductivity [ZZl I185J 
and surface roughness corrections |186| I187J . and obtained agreement with theory at 
better than 0.5% at the smallest separations of about 75 nm. However, potential 
corrections of greater than 1% exist, so that limits the level of verification of the 
theory. Their experimental work, which now continues at Harvard, suggests novel 
nanoelectromechanical applications. 

There is only one experiment with a parallel-plate geometry |188j . which is of 
limited accuracy (~ 15%) due to the difficulty of maintaining parallelism. It is, however, 
of considerable interest because the interpretation does not depend on the proximity 
theorem, corrections to which are problematic |148| l23j: see section l331 The importance 
of improving the accuracy of the parallel-plate configuration has been emphasized by 
Onofrio f^. 

The most precise experiment to date, using a MEMS, makes use of both static 
and dynamical procedures and yields a claimed accuracy of about 0.25% |;10m I190J . 
but this accuracy has been disputed jl91j . due to difficulty in controlhng roughness 
and the concomittant uncertainty in the ability to determine the separation distance. 
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It has been asserted |19(Jj that this experiment rules out the temperature dependence 
claimed in Ref. |lU9j (see section 13. 4.2j) . but this is problematic at this point, especially 
as comparison is only made with the MIM model ()H.53j) . rather than with the detailed 
calculation given there. 

Very recently, the Harvard group has performed a very interesting Casimir force 
measurement between a gold-covered plate and a sphere coated with a hydrogen- 
switchable mirror |192j . Although the mirror becomes transparent in the visible upon 
hydrogenation, no effect was observed on the Casimir force when the mirror was switched 
on and off. This shows that, in contradiction to the claims of Mostepanenko et al 
, for example in Ref. |127j . the Casimir force is responsive to a very wide range of 
frequencies, in accordance with the Lifshitz formula and the general dispersion relation 
for the permittivity. tl (See also Ref. |198| 1194] .) In particular, their results show that 
wavelengths much larger than the separation between the surfaces play a crucial role. 

Because all the recent experiments measure forces between relatively thin films, 
rather than between bulk metals, significant deviations from the Lifshitz formula (~ 2%) 
may be expected |195j . This may also be relevant to the claimed accuracy of the first 
Mohideen experiment [HZj, which uses a thin metallic coating, regarded as completely 
transparent. 

This may be an appropriate point to comment on the recent paper of Chen et al 
|95j . This is based on a reanalysis of experimental data obtained four years ago in 
Ref. jnni- Experimental precision of 1.75% and theoretical accuracy of 1.69% is claimed 
at the shortest distances, 62 nm. However, their analysis seems flawed. They obtain 
average experimental forces by averaging many measurements, which is only permissible 
if the averaging is carried out at exactly the same separation between the surfaces. Of 
course they have no way of knowing this. Furthermore, they apparently use the mean 
separation parameter rfo as a free variable in their fit, which essentially negates the 
possibility of testing the theory, which is most sensitive at the shortest separations |18()j . 
lannuzzi asserts that at distance of order 100 nm, errors of a few Angstroms preclude 
a 1% measurement. Therefore this analysis cannot be used as a serious constraint for 
either new forces or for setting limits on temperature corrections. ff 

A difference force experiment has been proposed by Mohideen and collaborators 
|196| I197J . The idea is to measure the difference in the force between a lens and a plate 
at room temperature, before and after both surfaces have been heated 50 K by a laser 
pulse. The measurements are not yet good enough to distinguish between the plasma 
and the Drude modes of the permittivity, or between the simplified impedance model 
versus the measured bulk permittivity approach, as discussed in section 13.41 

A proposal has been made to measure the force between eccentric cylinders, in which 
the axes are parallel but slightly offset |167j . The net force on the inner cylinder is zero, 

tt lannuzzi quotes Klinichitskaya as now agreeing with this statement. This, however, is hard to 
reconcile with statements made in Ref. |126| that one should use the extrapolated surface impedance 
value at lOc — l/2a rather than the actual zero-frequency value. 
ff I thank Davide lannuzzi for discussion of these points. 
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of course, when the cyhnders are concentric, but this equihbrium point is unstable. The 
idea is to look for a shift in the mechanical resonant frequency of the outer cylinder due 
to the Casimir force exerted by the inner one. The chief difficulty may be in maintaining 
parallelism. 

Another active area of experimental effort involving Casimir measurements is the 
search for new forces at the submicron level. These are based on looking for a discrepancy 
between the measured and predicted Casimir forces. The most recent limits are given 
in Krause, Decca, et al |198[ IKK?] . Unfortunately, the limits, for an assumed potential 
of the form 

\/(r) = -^^(l + ae-^/^), (3.92) 

for A < 10^^ m are only for absurdly large strengths, a < 10^^, and as A decreases 
the upper limit on a increases. The Purdue group has also proposed iso-electronic 
experiments to look at the force between a sphere and two different plates, composed 
of material with similar electronic properties (and hence similar Casimir forces) but 
different nuclear properties (and hence presumably different new forces). See Ref. |199j 
for a brief description of their experiment and the detection of a small, but probably 
not significant, residual force. 

Very recently, there has been a report |2()()j of an experiment |2()lj of dropping 
ultracold neutrons onto a surface. They are trapped between the mirror and the earth's 
gravitational potential. These gravitational bound states would be modified by any 
deviation from Newtonian gravity. No such deviations from Newton's law is found down 
to the 1-10 nm range. See also Nesvizhevsky and Protasov |2U2j who obtain limits on 
non-Newtonian forces inferior to those of Casimir measurements, that is, a < 10^^ at 
A = 10^^ m, although it is relatively better that the Casimir limits in the nanometer 
range, but the limits are extremely weak there, a < 10^^. 

It is clear that as micro engineering comes into its own, Casimir forces will have 
to be taken into account and utilized. A recent interesting paper by Chumak, Milonni, 
and Berman |2()3j suggests that the noncontact friction observed by Stipe et al |2()4j on 
a cantilever near a surface is due in major part to Casimir forces. The Casimir force 
is responsible for the frequency shift observed of about 4.5% for a gold sample at a 
separation of 2 nm. 

For another example along these lines, Lin et al |2n5j have shown that Casimir- 
Polder forces between atoms and the surface can provide fundamental limitations on 
stability of a Bose-Einstein condensate near a microfabricated silicon chip, a system 
which holds great promise for technological applications. 

The recent intense experimental activity is very encouraging to the development of 
the field. Coming years, therefore, promise ever increasing experimental input into a 
field that has been dominated by theory for five decades. 
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4. Self-Stress 

4-1- Surface and Volume Divergences 

It is well known that in general the Casimir energy density diverges in the neighborhood 
of a surface. For flat surfaces and conformal theories (such as the conformal scalar theory 
considered in Ref. jHZI, or electromagnetism) those divergences are not present.:}: We 
saw hints of this in section 12.41 In particular, Brown and Maclay |2n9j calculated the 
local stress tensor for two ideal plates separated by a distance a along the z axis, with 
the result for a conformal scalar 

^ ' 1440a4^ ^ ^ ^ ' 

This result was given recent rederivations in |2im Wf\ . Dowker and Kennedy |211j and 
Deutsch and Candelas |212j considered the local stress tensor between planes inclined 
at an angle a, with the result, in cylindrical coordinates (t, r, 6', z), 
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where for a conformal scalar, with Dirichlet boundary conditions, 

^<") = 2^ (^ - ^) • (") 

and for electromagnetism, with perfect conductor boundary conditions. 

For a — i> we recover the pressures and energies for parallel plates, (12.25 6j) . ()2.39j) and 
()3.87|) . (These results were later discussed in Ref. |213j .) 

Although for perfectly conducting flat surfaces, the energy density is flnite, for 
electromagnetism the individual electric and magnetic fields have divergent RMS values, 

(i^') ~ -(5') ~ ^, e^O, (4.5) 

a distance e above a conducting surface. However, if the surface is a dielectric, 
characterized by a plasma dispersion relation ()3.33|1 . these divergences are softened 

{E')r^\, -{B')r^\, e^O, (4.6) 



so that the energy density also diverges |214j 



(T°0) ~ i e ^ 0. (4.7) 

I In general, this need not be the case. For example, Romeo and Saharian j2(J6j show that with 
mixed boundary conditions the surface divergences need not vanish for parallel plates. For additional 
work on local effects with mixed (Robin) boundary conditions, applied to spheres and cylinders, and 
corresponding global effects, see Refs. [203 EOHl EOl • See also section El and Ref. [T^lTHj . 
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The null energy condition {rif^nf^ = 0) 

T^'n^n, > (4.8) 

is satisfied, so that gravity still focuses light. 

Graham |215j examined the general relativistic energy conditions required by 
causality. In the neighborhood of a smooth domain wall, given by a hyperbolic tangent, 
the energy is always negative at large enough distances. Thus the weak energy condition 
is violated, as is the null energy condition ()4.8p . However, when ()4.8p is integrated over 
a complete geodesic, positivity is satisfied. It is not clear if this last condition, the 
Averaged Null Energy Condition, is always obeyed in flat space. Certainly it is violated 
in curved space, but the effects always seem small, so that exotic effects such as time 
travel are prohibited. 

However, as Deutsch and Candelas |212j showed many years ago, in the 
neighborhood of a curved surface for conformally invariant theories, (T^;^) diverges as 
e~^, where e is the distance from the surface, with a coefficient proportional to the sum 
of the principal curvatures of the surface. In particular they obtain the result, in the 
vicinity of the surface, 

(r,.)~e-^r(3) + e-2T(2,) + e-^r«, (4.9) 

and obtain explicit expressions for the coefficient tensors T^J and T^j, in terms of the 
extrinsic curvature of the boundary. 

For example, for the case of a sphere, the leading surface divergence has the form, 
for conformal fields, f or r = a + e, e ^ 

f 2/a \ 
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in spherical polar coordinates, where the constant is A = l/14407r^ for a scalar, or 
A = l/1207r^ for the electromagnetic field. Note that ()4.10|) is properly traceless. The 
cubic divergence in the energy density near the surface translates into the quadratic 
divergence in the energy found for a conducting ball |216j . The corresponding quadratic 
divergence in the stress corresponds to the absence of the cubic divergence in (T^r) • 

This is all completely sensible. However, in their paper Deutsch and Candelas 
|212j expressed a certain skepticism about the validity of the result of Ref. [23 for 
the spherical shell case (described in part in section 14. 4|) where the divergences cancel. 
That skepticism was reinforced in a later paper by Candelas |217j . who criticized the 
authors of Ref. |21] for omitting 6 function terms, and constants in the energy. These 
objections seem utterly without merit. In a later critical paper by the same author |218j . 
it was asserted that errors were made, rather than a conscious removal of unphysical 
divergences. 



{T^Mu) - 3 



(4.10) 
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Of course, surface curvature divergences are present. As Candelas noted |217tl218j . 
they have the form 

E = E^ f dS + E^ f dS {ki + K2) + Ef I dS {ki - K^f + E^j f dSKiK2 + ..., (4.11) 

where Ki and K2 are the principal curvatures of the surface. The question is to what 
extent are they observable. After all, as has been shown in Ref. [211 EZI and in section|231 
we can drastically change the local structure of the vacuum expectation value of the 
energy-momentum tensor in the neighborhood of flat plates by merely exploiting the 
ambiguity in the definition of that tensor, yet each yields the same finite, observable 
(and observed!) energy of interaction between the plates. For curved boundaries, much 
the same is true. A priori, we do not know which energy-momentum tensor to employ, 
and the local vacuum-fluctuation energy density is to a large extent meaningless. It 
is the global energy, or the force between distinct bodies, that has an unambiguous 
value. It is the belief of the author that divergences in the energy which go like a power 
of the cutoff are probably unobservable, being subsumed in the properties of matter. 
Moreover, the coefficients of the divergent terms depend on the regularization scheme. 
Logarithmic divergences, of course, are of another class [21] • 

Dramatic cancellations of these curvature terms can occur. It might be thought 
that the reason a finite result was found for the Casimir energy of a perfectly conducting 
spherical shell |19| |221 121] is that the term involving the squared difference of curvatures 
in ()4.11|) is zero only in that case. However, for reasons not yet apparent to the 
present author, it has been shown that at least for the case of electromagnetism the 
corresponding term is not present (or has a vanishing coefficient) for an arbitrary smooth 
cavity |219j , and so the Casimir energy for a perfectly conducting ellipsoid of revolution, 
for example, is finite. This finiteness of the Casimir energy (usually referred to as the 
vanishing of the second heat-kernel coefficient ^U]) for an ideal smooth closed surface 
was anticipated already in Ref. [22], but contradicted by Ref. |212j . More specifically, 
although odd curvature terms cancel inside and outside for any thin shell, it would be 
anticipated that the squared-curvature term, which is present as a surface divergence in 
the energy density, would be refiected as an unremovable divergence in the energy. 
For a closed surface the last term in ()4.1H) is a topological invariant, so gives an 
irrelevant constant, while no term of the type of the penultimate term can appear 
due to the structure of the traced cylinder expansion [70] . It would be extraordinarily 
interesting if this Casimir energy could be computed for an ellipsoidal boundary, but the 
calculation appears extremely difficult because the Helmholtz equation is not separable 
in the exterior region. 

4-2. Casimir Forces on Spheres via S- Function Potentials 

This section is an adaptation and an extension of calculations presented in Ref. |62j . 
This investigation was carried out in response to the program of the MIT group 
[165| Wi \ l?)T | I7T | IH( )j . They rediscovered irremovable divergences in the Casimir energy for 



The Casimir Effect 



41 



a circle in 2 + 1 dimensions first discovered by Sen f2'20\ I221J , but then found divergences 
in tlie case of a splierical surface, thereby casting doubt on the vahdity of the Boyer 
calculation (TH]. Some of their results, as we shall see, are spurious, and the rest are 
well known ^^. However, their work has been valuable in sparking new investigations 
of the problems of surface energies and divergences. 

We now carry out the calculation we presented in section |21 in three spatial 
dimensions, with a radially symmetric background 
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which would correspond to a Dirichlet shell in the limit A 
transformed Green's function satisfies the equation (k^ 
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We write Q in terms of a reduced Green's function 
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oo. The time- Fourier 

(4.13) 
(4.14) 
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gi{r,r') = —S{r-r'). (4.15) 
r^ a ' ' 

We solve this in terms of modified Bessel functions, Iu{x), K,y{x), where u = I + 1/2, 
which satisfy the Wronskian condition 
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The solution to ()4.15|) is obtained by requiring continuity of gi at each singularity, r' 
and a, and the appropriate discontinuity of the derivative. Inside the sphere we then 
find (0 < r,r' < a) 
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Here we have introduced the modified Riccati-Bessel functions, 

si{x) = J—Ii+i/2{x), ei{x) = J—Ki+i/2{x). 



(4.17) 
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Note that ()4.17p reduces to the expected result, vanishing as r — ^ a, in the limit of 
strong coupling: 
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When both points are outside the sphere, r,r' > a, we obtain a similar result: 
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which similarly reduces to the expected result as A — > oo. 
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Now we want to get the radial-radial component of the stress tensor to extract the 
pressure on the sphere, which is obtained by applying the operator 



drdr>-^i-d'd" + V-V')^^ 



drdr' — K^ — 






(4.21) 



to the Green's function, where in the last term we have averaged over the surface of the 
sphere. In this way we find, from the discontinuity of (Trr) across the r = a surface, the 
net stress 

1=0 -^^ ^ ^ X 

The same result can be deduced by computing the total energy (|2.15p . The free 
Green's function, the first term in ()4.17p or ()4.20j) . evidently makes no significant 
contribution to the energy, for it gives a term independent of the radius of the sphere, 
a, so we omit it. The remaining radial integrals are simply 

j dye}{y) = -—[{x^ + l{l + l))e^i+xeie[-x\f], (4.236) 

where all the Bessel functions on the right-hand-sides of these equations are evaluated 
at X. Then using the Wronskian, we find that the Casimir energy is 

1 °° /"°° d 

E = -^—Y{2l + l) dxx — \Yi[l + \h{x)K,{x)]. (4.24) 

27ra -^ Jo dx 

If we differentiate with respect to a, with A/a fixed, we immediately recover the force 
()4.22|) . This expression, upon integration by parts, coincides with that given by Barton 
|222j , and was first analyzed in detail by Scandurra |223j . For strong coupling, it reduces 
to the well-known expression for the Casimir energy of a massless scalar field inside and 
outside a sphere upon which Dirichlet boundary conditions are imposed, that is, that 
the field must vanish at r = a: 

1 °° r^ d 

lim E = -^^y(2/ + l) / Axx — \n[I^{x)K^{x)], (4.25) 

because multiplying the argument of the logarithm by a power of x is without effect, 
corresponding to a contact term. Details of the evaluation of Eq. fj4.25|) are given in 
Ref. IHZI, and will be considered in section |13l below. (See also Refs. p2^ l225l I226J .1 

The opposite limit is of interest here. The expansion of the logarithm is immediate 
for small A. The first term, of order A, is evidently divergent, but irrelevant, since that 
may be removed by renormalization of the tadpole graph. In contradistinction to the 
claim of Refs. jHH IHOl IHOl El, the order A^ term is finite, as established in Ref. [FTj . 
That term is 

^^''^ = i^E(2^ + l)/ dxx-^[h^,f,{x)K,^y,{x)]'. (4.26) 

^ 1=0 "^0 



The Casimir Effect 43 

The sum on / can be carried out using a trick due to Klich |227j : The sum rule 

oo 

V(2/ + l)ei{x)si{y)Pi{cose) = ^e"", (4.27) 



where p = \/ x^ + y^ — 2xy cos 6*, is squared, and then integrated over 6*, according to 

I d(cos^)PKcos^)P,.(cos^)=(5«,^^. (4.28) 

In this way we learn that 

5^(2/ + l)e,2(x).f(x) = ^/ ^e-. (4.29) 

1=0 '"^ 

Although this integral is divergent, because we did not integrate by parts in ()4.26|) . that 
divergence does not contribute: 

E^^> = - — / dx-x-- — e "" = , (4.30) 

Ana Jq 2 dx Jo "^ 327ra 

which is exactly the result (4.25) of Ref. jSZ], which also follows from ()2.2()p here. 

However, before we wax too euphoric, we recognize that the order A^ term appears 

logarithmically divergent, just as Refs. jHOI and ^J claim. This does not signal a 

breakdown in perturbation theory, as the divergence ()2.2H) in the D = 1 calculation 

did. Suppose we subtract off the two leading terms. 
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To study the behavior of the sum for large values of /, we can use the uniform asymptotic 
expansion (Debye expansion). 
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Here x = uz, and t = 1/Vl + z^- The functions A and B, etc., are polynomials in t. 
We now insert this into ()4.3H) and expand not in A but in u; the leading term is 

^'^''-^Eif(r^-^c(i). (4.33) 

Although the frequency integral is finite, the angular momentum sum is divergent. The 
appearance here of the divergent ({1) seems to signal an insuperable barrier to extraction 
of a finite Casimir energy for finite A. The situation is different in the limit A -^ oo - 
See section lOl 

This divergence has been known for many years, and was first calculated explicitly 
in 1998 by Bordag et al [SHI, where the second heat kernel coefficient gave an equivalent 
result, 

\3 ^ 

E , s^O. (4.34) 
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A possible way of dealing with this divergence was advocated in Ref. |223j . Very 
recently, Bordag and Vassilevich |228j have reanalyzed such problems from the heat 
kernel approach. They show that this O(A^) divergence corresponds to a surface tension 
counterterm, an idea proposed by me in 1980 [23 I229J in connection with the zero- 
point energy contribution to the bag model. Such a surface term corresponds to A/a 
fixed, which then necessarily implies a divergence of order A^. Bordag argues that it 
is perfectly appropriate to insert a surface tension counterterm so that this divergence 
may be rendered finite by renormalization. 



4-3. TM Spherical Potential 

Of course, the scalar model considered in the previous subsection is merely a toy model, 
and something analogous to electrodynamics is of far more physical relevance. There are 
good reasons for believing that cancellations occur in general between TE (Dirichlet) and 
TM (Robin) modes. Certainly they do occur in the classic Boyer energy of a perfectly 
conducting spherical shell f^l 1221211! and the indications are that such cancellations 
occur even with imperfect boundary conditions |222j . Following the latter reference, let 
us consider the potential 
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In the limit A — > oo this corresponds to TM boundary conditions. The reduced Green's 
function is thus taken to satisfy 
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At r = r' we have the usual boundary conditions, that gi be continuous, but that its 
derivative be discontinuous, 



2d 
or 

while at the surface of the sphere the derivative is continuous, 
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while the function is discontinuous, 
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(4.38a) 
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Equations ()4.38a|) and ()4.386|) are the analogues of the boundary conditions ()2.30a|) . 
fl2.306|) treated in section [231 

It is then easy to find the Green's function. When both points are inside the sphere. 
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and when both points are outside the sphere. 
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It is immediate that these supply the appropriate Robin boundary conditions in the 
A ^ CX3 hmit: 
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The Casimir energy may be readily obtained from ()2.15|) . and we find, using the 
integrals (g^Sil), (14123^ 

1 °° f°° d 

E = V(2Z + 1) / dxx—\n[l + \xe'(x)s'(x)]. (4.41) 
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The stress may be obtained from this by applying —d/da, and regarding Xa as constant 
[see ()4.35|1 ]. or directly, from the Green's function by applying the operator. 
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which is the same as that in ()4.2ip . except that 



-drV, 



(4.42) 



(4.43) 



appropriate to TM boundary conditions (see Ref. |23()j . for example). Either way, the 
total stress on the sphere is 
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27ra2 -^ ./n 1 + Xxe'i{x) s[{x) 

The result for the energy (j4.4ip is similar, but not identical, to that given by Barton 



Suppose we now combine the TE and TM Casimir energies, ()4.24|1 and ()4.41|1 : 
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In the limit A — »• cxo this reduces to the familiar expression for the perfectly conducting 
spherical shell |24j : 



lim E = — y (2/ + 1)/ dxx(-^ + % + ^ + -\\. (4.46) 

A^oo 27ra -^ Jo \^i e, si sj 

Here we have, as appropriate to the electrodynamic situation, omitted the I = mode. 
This expression yields a finite Casimir energy, as we will see in section I^Til What about 
finite A? In general, it appears that there is no chance that the divergence found in the 
previous section in order A^ can be cancelled. But suppose the coupling for the TE and 
TM modes are different. If A^^A"""^ = 4, a cancellation appears possible. 

Let us illustrate this by retaining only the leading terms in the uniform asymptotic 
expansions: (x = uz) 
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Then the logarithm appearing in the integral for the energy (j4.45|) is approximately 

m r. In ^-^ +^-[^^^y 1- (^ - A^j • (^•^«) 

The first term here presumably gives no contribution to the energy, because it is 
independent of A upon differentiation, and further we may interpret Yl'i^o ^^ = ^ [see 
()4.52|) ]. Now if we make the above identification of the couplings, 

\TE 2 

^ = ^ = Xtm' ^4.49) 

all the odd powers of z/ cancel out, and 

1-^1. (4.50) 
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E^ V(2/ + l)/ dxx— In 

27ra ^ Jo dx 



v^ 



The divergence encountered for the TE mode is thus removed, and the power series is 
simply twice the sum of the even terms there. This will be finite. Presumably, the same 
is true if the subleading terms in the uniform asymptotic expansion are retained. 

It is interesting to approximately evaluate (J4.5(J|) . The integral over z may be easily 
evaluated as a contour integral, leaving 

This / sum appears to be divergent, an artifact of the asymptotic expansion, since we 
know the A^ term is finite. However, if we expand the square root for small A^/z/^, we 
see that the O(A^) term vanishes if we interpret the sum as 

oo 

5^z.- = (2^-l)C(.), (4.52) 

/=o 

in terms of the Riemann zeta function. The leading term is O(A^): 

A^ ^ 1 AV2 

E > ^ = . 4.53 

z=o 
To recover the correct leading A behavior in ()4.3U|) requires the inclusion of the 
subleading z/"^" terms displayed in ()4.32p . 

Much faster convergence is achieved if we consider the results with the / = term 
removed, as appropriate for electromagnetic modes. Let's illustrate this for the order 
A^ TE mode (now, for simplicity, write A = A"""^) Then, in place of the energy (j4.3Up we 
have 

E(^^) = ^^^r% sinh^ xe-^^ = - f — + ^1 = -(0.065 106 1). (4.54) 

Now the leading term in the uniform asymptotic expansion is no longer zero: 



EW = -^y(2/+i) / dxx^ 

2'Ka ^-^ In ax 

1=1 ^^ 

1,2 °o ^ ^x \2 \2 



^Z/2 



5:^°(-^)-1^->-0625). (4.55) 



Svra ^-^ \ 2 / 16a a 
1=1 
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which is 4% lower than the exact answer (J4.54J) . The next term in the uniform asymptotic 
expansion is 



E(2) = __^[3^(2)_4] / dzf 
Ana' Vo 



8 



= ^ (^ - — 1 = -(0.002 736 8), (4.56) 

a V2048 256y a ^ ^' ^ ^ 

which reduces the estimate to 

^(0) _^ ^(2) _ _(o.065 236 8), (4.57) 

a 

which is now 0.2% high. Going out one more term gives 

\2 roo J.4 

^(4) = _ J^ [15^(4) _ 16] / dzf — {7 - U8t^ + 554t^ - 708t^ + 295t«) 

•-' u 

A2 / 597r4 177 \ A^ , , , 

= = 0.000 158 570 , 4.58 

a V524288 16328/ a ^ ^' ^ ^ 

and the estimate for the energy is now only 0.04% low: 

\2 

^(0) ^ ^(2) ^ ^(4) ^ _(o.065 078 23). (4.59) 

a 

We could also make similar remarks about the TM contributions. 

4.4- Perfectly Conducting Spherical Shell 

Now we consider a massless scalar in three space dimensions, with a spherical boundary 
on which the field vanishes. This corresponds to the TE modes for the electrodynamic 
situation first solved by Boyer [13 E21 121|- The purpose of this section (adapted from 
Ref. |57jj) is to emphasize anew that, contrary to the imphcation of Ref. [UUIIHH 1301 171] . 
the corresponding Casimir energy is also finite for this configuration. 

The general calculation in D spatial dimensions was given in Ref. |224j ; the pressure 
is given by the formula 

p-t ',';:;.;. Jrl'^iK-- f'^'^^ '■' [^"<^'^^"(^)^^"''] ■ (^■»' 

Here z/ = / — 1 + D/2. For D = 3 this expression reduces to 

1 °° f°^ d 

^ = -^^^^(2^ + 1)/ dxx—\n[li+y2ix)Ki+y2ix)/x]. (4.61) 

This precisely corresponds to the strong limit A — > oo given in ()4.25j) . if we recall the 
comment made about contact terms there. In Ref. |224j we evaluated expression ()4.60p 
by continuing in D from a region where both the sum and integrals existed. In that 
way, a completely finite result was found for all positive D not equal to an even integer. 
Here we will adopt a perhaps more physical approach, that of allowing the time- 
coordinates in the underlying Green's function to approach each other, as described in 
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Ref. p4j. That is, we recognize that the x integration above is actually a (dimensionless) 
imaginary frequency integral, and therefore we should replace 

/ dxfix) = - dye'y'fi\y\), (4.62) 

Jo ^ J~oo 

where at the end we are to take 5 — *> 0. Immediately, we can replace the x~^ inside the 
logarithm in ()4.61|) by x, which makes the integrals converge, because the difference is 
proportional to a delta function in the time separation, a contact term without physical 
significance. 

To proceed, we use the uniform asymptotic expansions for the modified Bessel 
functions, ()4.32j) . This is an expansion in inverse powers of z/ = / + 1/2, low terms in 
which turn out to be remarkably accurate even for modest /. The leading terms in this 
expansion are 

In [xIi+i/2{x)Ki+y2{x)\ ~ In y + —g{t) + —h(t) + . . . , (4.63) 

where x = vz and t = (1 + 2^)~^/^. Here 

g{t) = ^-{f-Qt' + ht'), (4.64a) 

h{t) = 7T(13t^ - 284t^ + 1062t® - 1356t^° + 565^^^). (4.64&) 

The leading term in the pressure is therefore 

1 °° poo -| °° q 

Po = -^^,J2(^l + l> dzt2 = --^5^z.2 = -^C(-2) = 0. (4.65) 

Svr^a* -^^ Jn Svra* -^^ 327ra* 

1=0 ■^'^ 1=0 

where in the last step we have used the formal zeta function evaluation (j4.52|) .^ Here 
the rigorous way to argue is to recall the presence of the point-splitting factor e"^^*^ and 
to carry out the sum on / using 

OO -, -I 



^^ 2ismzS/2' ^ ' 



1=0 

so 

Vz/V^^^ = - ^ = -( ^ + i ^ (4 67) 

^ d{z5Y2siYiz5/2 S\ sm^z5/2 smz6/2j' ^' ' 

Then Pq is given by the divergent expression 

which we argue is zero because the integrand is odd, as justified by averaging over 
contours passing above and below the pole at 2; = 0. 

The next term in the uniform asymptotic expansion (j4.63p . that involving g, likewise 
gives zero pressure, as intimated by the formal zeta function identity ()4.52j) . which 

§ Note that the corresponding TE contribution for the electromagnetic Casimir pressure would not be 
zero, for there the sum starts from / = 1. 
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vanishes at s = 0. The same conclusion follows from point splitting, using (J4.66J) and 
arguing that the resulting integrand ~ z^t^ g' {t) / z5 is odd in z. Again, this cancellation 
does not occur in the electromagnetic case because there the sum starts at / = 1. 
So here the leading term which survives is that of order v~^ in ()4.63p . namely 



CO _ 






dzh{t), 



(4.69) 



where we have now dropped the point-splitting factor because this expression is 
completely convergent. The integral over z is 

357r 



dzhit) 



(4.70) 



,0 32768 

and the sum over / is 3^(2) = 7r^/2, so the leading contribution to the stress on the 
sphere is 

357r2 0.00527094 



^2 = 47ra^P2 



65536a2 



(4.71) 



Numerically this is a terrible approximation. 

What we must do now is return to the full expression and add and subtract the 
leading asymptotic terms. This gives 



^ CO 



1=0 



where 



Ri = Qi + 



dx 



\nzt + \g{t) + \h{t) 



V 



where the integral 

Qi = - 



dx \v^xIv{x^Ky{x)\ 



was given the asymptotic form in Ref. 
z/vr vr 357r 5657r 



Q 



2 ^ 128z/ 



(/ » 1): 
12087677r 



1380083577r 



(4.72) 



(4.73) 



(4.74) 



(4.75) 



32768Z/3 1048577Z/5 2147483648//^ 137438953472z/9 ' 
The first two terms in ()4.75p cancel the second and third terms in ()4.73|) . of course. 
The third term in (J4.75J) corresponds to /i(t), so the last three terms displayed in (|4.75j) 
give the asymptotic behavior of the remainder, which we call wiv). Then we have, 
approximately. 



^ n ^ oo 

S^S^ {y^vRi ^ V vwiv). (4.76) 

/=o l=n+\ 

For n = \ this gives S ~ 0.002 852 78/a^, and for larger n this rapidly approaches the 
value first given in Ref. |224j . and rederived in |225t 12261 231] 

5™ = 0.002817/a^ (4.77) 
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a value much smaller than the famous electromagnetic result fHl 1^ 1^ 122] , 



0.04618 



S^^ = 7 , (4.78) 



a 



because of the cancellation of the leading terms noted above. Indeed, the TM 
contribution was calculated separately in Ref. |230j . with the result 

S™ = -0.02204^, (4.79) 



a 



and then subtracting the / = modes from both contributions we obtain ()4.78|) 

48a2 a^ 

4-5. Dielectric Spheres 

The Casimir self-stress on a uniform dielectric sphere was first worked out in 1979 
|216j . It was generalized to the case when both electric permittivity and magnetic 
permeability are present in 1997 |2H2j . Since this calculation is summarized in my 
monograph ^25^^ , we content ourselves here with simply stating the result for the pressure, 
(x = y/efl\y\, x' = y/e'fj,'\y\ where e',fi' are the interior, and e,fi are the exterior, values 
of the permittivity and the permeability) 

1 r dy .. ^ 2/ + 1 r d , ^ 
p = — -/ —e'y^y <^x— InA 



2a^ y oo 2vr ^^ 47r 1 dx 
+ 2x'[s\ix')eKx') - e,ix')s'lix')] - 2x[s[ix)e\ix) - Q(x).r(x)]| , (4.81) 

where the "bulk" pressure has been subtracted, and 

A = [si{x')e[ix) - s[{x')e,{x)r - e[siix')e[{x) + s[{x')ei{x)]\ (4.82) 

with the parameter ^ being 



e' ^l 



e = ^^ , (4.83) 

7^+1 

and S is the temporal regulator introduced in ()4.62|) . This result is obtained either by 
computing the radial-radial component of the stress tensor, or from the total energy. 



In general, this result is divergent. However, consider the special case y^ = -^/e'/i', 
that is, when the speed of light is the same in both media. Then x = x' and the Casimir 
energy derived from (j4.8ip reduces to 

E = Ana^P = "^ / dye'y' ^(2/ + l)x— ln[l - eUsieiYf], (4.84) 



1=1 
where 

_ n- jj,' _ e-e' 
fi + fi' e + e' 



(4.85) 
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If ^ = 1 we recover the case of a perfectly conducting spherical shell, treated in section lOl 
[cf. flOHll ]. for which E is finite. In fact fHl^ is finite for all ^. 
Of particular interest is the dilute limit, where |227j 
5^2 0.099 4718^2 

[This evaluation is carried out in the same manner as that of ()4.26|) .] It is remarkable 
that the value for a spherical conducting shell ()4.78p . for which ^ = 1, is only 7% lower, 
which as Klich remarks, is accounted for nearly entirely by the next term in the small 
^ expansion. 

There is another dilute limit which is also quite surprising. For a purely dielectric 
sphere (/^ = 1) the leading term in an expansion in powers of e — 1 is finite 
p33l 123^ 15911235]: 

^.^(£^.(,,l)^^£21Z5!. (4.87) 

ISdovr a a 

This result coincides with the sum of van der Waals energies of the material making up 

the ball ^^. The term of order (e — 1)^ is divergent [221 • The establishment of the result 

(I4.87P was the death knell for the Casimir energy explanation of sonoluminescence |23(jj 

- See section ini 

The temperature correction to this result was first worked out by Nesterenko, 

Lambiase, and Scapetta |237| 1238] . See also Ref . |239j . 

^.^. Cylinders 

It is much more difficult to carry out Casimir calculations for cylindrical geometries. 
We restrict our attention here to cylinders of circular cross section and infinite 
length. Although calculations have been carried out for parallelopiped geometries 
pm Ym\. Yim Yim Yim. Yim Yim. YWA Yim Yim l^nni YIF^ . the effects included refer 
only to the interior modes of oscillation. This is because the wave equation is not 
separable outside a cube or a rectangular solid. As a result, divergences occur which 
cannot be legitimately removed, which nevertheless are artificially removed by zeta- 
function methods. It is the view of the author that such finite results are without 
meaning. 

But even though circular-cylinder calculations are possible, they are considerably 
more complex than the corresponding spherical calculations. This is not merely 
because spherical Bessel functions are simpler than cylinder functions. The fundamental 
difficulty in these geometries is that there is in general no decoupling between TE and 
TM modes j252j . Progress in understanding has therefore been much slower in this 
regime. It was only in 1981 that it was found that the electromagnetic Casimir energy 
of a perfectly conducting cylinder was attractive, the energy per unit length being j253] 



^em,cyl — 5 , (4.88) 



a? 
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for a circular cylinder of radius a. The corresponding result for a scalar field satisfying 
Dirichlet boundary conditions of the cylinder is repulsive 



0.000606 ,,„„, 

^D,cyi = 5 • 4.89 

These ideal limits are finite, but, as with the spherical geometry, less ideal 
configurations have unremovable divergences. For example, a cylindrical 5-shell 
potential, as described earlier, has divergences (in third order) |255j . And it is expected 
that a dielectric cylinder will have a divergent Casimir energy, although the coefficient 
of (e — 1)^ will be finite for a dilute dielectric cylinder |256j . corresponding to a finite 
van der Waals energy between the molecules that make up the material |257j . Recent 
progress in understanding these points will be described below. 

4-6.1. Dielectric cylinders The following calculation represents work in progress with 
Ines Cavero-Pelaez. Although the calculation remains incomplete, we offer it here as 
a detailed example of how a complicated electromagnetic calculation is formulated in 
the Green's function approach. We start from the equations satisfied by the Green's 
dyadics for Maxwell's equations in a medium characterized by a permittivity e and a 
permeability /i (see Ref. |216j ): 

1 
V X r' - iujfi^ = -V X 1, (4.90a) 



-V X $ - itueT' = 0, (4.906) 



where 



r'(r, r', Lu) = r(r, r'; to) + ^- (4.91) 

and where the unit dyadic 1 includes a three-dimensional 6 function, 

l = lS{r-r'). (4.92) 

The two dyadics are solenoidal, 

V ■ $ = 0, (4.93a) 

V-r' = 0. (4.936) 

The corresponding second-order equations are 

(V^ + uj'^efi)T' = --V X (V X 1), (4.94a) 

£ 

(V^ + cj V)^ = i^V X 1. (4.946) 

Classically, these Green's dyadic equations are equivalent to Maxwell's equations, 
and give the solution thereto when a polarization source P is present, 

^W^/(d.')r(.,.')-PM. (4.95, 
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Quantum mechanically, they give the one-loop vacuum expectation values of the product 
of fields (at a given frequency uj) 

{E{r)E{r')) =-T{r,r'), (4.96a) 

(H(r)H(r')) = --^-V x T(r, r') x V. (4.966) 

1 UJ^fi^ 

Thus, from knowledge of the classical Green's dyadics, we can calculate the one-loop 
vacuum energy or stress. 

We now introduce the appropriate partial wave decomposition for a cylinder, a 
slight modification of that given for a conducting cylindrical shell 



'~j<J /*00 1 7 

T'{r,r';uj)= ^ / —l^{Vxz)f^{r;k,uj)xmk{0,z) 

+ — Vx {V xz)g^{r;k,iu)xmkiO,z)}, (4.97a) 

°2^ /"°° dA; f 

^{r,r';uj)= 2_, / ^[{^ ^ z)gm{r;k,uj)xmk{0,z) 

—V X (V X £)/^(r; k, iu)xmki0, z)], (4.976) 

UJU J 



LUfI 

where the cylindrical harmonics are 

1 



Xmk{0,z) = ^=e'^'e''\ (4.98) 

V Svr 



dm = z—J— - — • (4-100) 



and the dependence of fm etc. on r' is implicit (they are further vectors in the second 
tensor index). Because of the presence of these harmonics, we have 

vm '^ o 

V X i ^ f 6— = M, (4.99a) 

r or 

V X (V X i) ^ fik— - 0— - zdra = Af, (4.996) 

or r 

in terms of the cylinder operator 
1 d d m? 

Now if we use the Maxwell equation ()4.906|) we conclude^ 

gm = gm, (4.101a) 

(rf™ - k^)fm = -UJ^l^fm. (4.1016) 

From the other Maxwell equation ()4.90a|l we deduce (we now make the second, 
previously suppressed, position arguments explicit) (the prime on the differential 

II It might be thought that we could immediately use the general waveguide decomposition of modes 

into those of TE and TM type, for example as given in Ref. j258| . However, this is here impossible 

because the TE and TM modes do not separate. See Ref. (252j . 

% The ambiguity in solving for these equations is absorbed in the definition of subsequent constants of 

integration. 
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operator signifies action on the second, primed argument) 

,,.2 

dm'Dm9m{r;r',e',z') 
where the Bessel operator appears, 

T^m = dm + A^, A^ = Lu'^EfJ. — k"^. 

Now we do the separation of variables on the second argument. 



e r 



U{r, r') = M*'FUr, r'; k, u) + J\f*'FUr, r'; k, u) 



Xmk{(^,z), 



gm[r,r 



1 

UJ 
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(4.102a) 
(4.1026) 

(4.103) 

(4.104a) 
(4.1046) 



N-*'GUr, r- k, Lu) + M*'GUr, r'- k, lu)\ x*mkiO', z'), 
where we have now introduced the two scalar Green's functions Fm, Gm, which satisfy 

CJ^/i 1 



(^m^m-'^myj ^ ) 



--6{r — r') 



e r 
2I 



dmVmGmir, v) = UJ -6{r - r'). 



r 



(4.105a) 
(4.1056) 



while Fm and Gm are annihilated by the operator dmT>m- 

In the following we will apply these equations to a dielectric-diamagnetic cylinder 
of radius a, where the interior of the cylinder is characterized by a permittivity e and 
a permeability /i, while the outside is vacuum, so e = fi = 1 there. Let us compute the 
Green's dyadics for the case that the source point is outside, r' > a. If the field point is 
also outside, r,r' > a, the Green's dyadics have the form {fi = e = 1) 



r' 



E fiW 



u^ 



M*'Fm{r, r'- k, to) + Ar*'Fm{r, r'; k, cu) 



1 



$ 



+ —AT Ar*'Gmir, r'- k, Lu) + M*'Gmir, r'- k, lu) \xmkiO, z)x*mkiO', z'), (4.106a 

00 



dA; r i 



27r leu 



M 



Af*'Gm{r, t'- k, u) + M*'Gm{r, r'; k, u) 



UJ 



-J\r M*'Fm{r, r'; k, u) + M*'Fm{r, r'; fc, u) 



Xmk{0,z)Xmki^',z'). 



(4.1066) 



From the differential equation ()4.1()5a| we see that the Green's function F has the form 

2 



F =-- 
A2 



2\m\ \r 



— ] + 7^Jm{kr<^)Hm{kr>) 



2i 



+ amHm{\r)Hm{\r') + 6„r-l'^li/„(Ar') + c„r'-l'"liy„(Ar) + t/„ 



'm\ i-\m\ 



,r ' T 



(4.107) 



while Gm has the same form with the constants am, bm, Cm, dm replaced by a^, 6J„, c^, 
d'j^, respectively. The homogeneous functions have the form 



Fm = dmHm{Xr)Hm{Xr') + hmr~^""^Hm{\r') + g™r'-l'"l/7„(Ar) + d„ 



\m\ t-\m\ 



r ' T 



, (4.108) 
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and Gm replaces a -^ a', etc. 

When the source point is outside and the field point is inside, there are only 
homogeneous solutions of the equations, so we may write for r < a,r' > a 

Fm = ey^\r'-\'^\ + fy'^\H„,{\r') + ^„ J„(A'r)r'~l™l + /i„ J„,(AV)iJ„(Ar'), (4.109) 

and similarly for G^, -^m, G^, with the constants denoted by e^, Cm, and e^, 
respectively. Here the outside and inside forms of A are given by 

X^=uo^-k\ \'^ = cu''fxe-k^ (4.110) 

The various constants are to be determined, as far as possible, by the boundary 
conditions at r = a. The boundary conditions at the surface of the dielectric cylinder are 
the continuity of the tangential components of the electric field, of the normal component 
of the electric displacement, of the normal component of the magnetic induction, and 
of the tangential components of the magnetic field: 

Et is continuous, eEn is continuous, 

Hf is continuous, fiHn is continuous. (4.111) 

These conditions are redundant, but we will impose all of them as a check of consistency. 
In terms of the Green's dyadics, the conditions read 

(4.112a) 
(4.112&) 
(4.112c) 
(A.lUd) 
(4.112e) 

(4.112/) 

A fairly elaborate system of linear equations for the various constants results. However, 
they are not quite sufficient to determine all the relevant physical combinations. We 
also need to impose one of the Helmholtz equations, say ()4.946|) . From that equation 
we learn 

b' — ksgnmb = 0, (4.113a) 

b-'-^b' = 0, (4.113?.) 

k 

d + d' = 0, (4.113 c) 

f^t?i^f' = 0, (4.113d) 

e k 

f' + -ksgnmf = 0, (4.113e) 

/i 

e-^e' = 0, (4.113/) 



e 


r 


T=a 


is continuous. 


z 


v 


r=a 


is continuous. 


r 


eV 


t 
r=a 


is continuous 


r 


^^^ 


> 

r=a 


is continuous 





$ 


r=a 


is continuous. 


z 


$ 


r=a 


is continuous. 
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where we have introduced the abbreviations for any constant K 
K = K -k sgnm K, K' = K' 



sgnm^, 



k 



56 



(4.114) 



Then from the boundary conditions we can solve for the remaining constants: First, 



and 



c = 
9 = 

h! 



^m 



0, 
^0, 



/i 



6fx) 



uj'^mk 



XX'D 

uj'^mk 



hmHm{Xa) Jm{X tt), 
hmJm{X a) , 






+ 



hr 



all in terms of 



X^2i H^{\a) A2/i '"iJ„(Aa) 



e S 



where the denominators occurring here are 

D = eXaJ'^{X'a)H^{Xa) - X'aJ^{X'a)H'^{Xa), 
D = fiXaJ'^{X'a)Hm{Xa) — X'aJm{X'a)H'^{Xa), 
E = {XX'fDD - {sfi - lfk^mWHl{Xa)Jl{X'a) 

The second set of constants is 



- — (1 - en)-—^h'^HmiXa)JmiX'a), 



XX'D 

A'^ (1 — Efi) mk 



A2 



AA'Z^ 



=:hj^Jm{X a) 



in terms of 



uj"^ -n JmjXa) A^ JmjX'a) 
X^2i H^{Xa) A% ""HrniXa) 



AX'D 



It might be thought that m = is a special case, and indeed 






lln^, 
2 r^ 



4.115a) 
4.1156) 

4.116a) 
4.1166) 
4.116 c) 

4:.116d) 



4.117a) 


4.1176) 


4.117c) 


4.118a) 


4.1186) 


4.118c) 



(4.118(i) 



(4.119) 



2\m\ \i>/ ^ /> 

but just as the latter is correctly interpreted as the limit as \m\ -^ 0, so the coefficients 
in the Green's functions turn out to be just the m = limits of those given above, so 
the 771 = case is properly incorporated. 

It is now easy to check that, as a result of the conditions ()4.113aj) . ()4.1136|) . ()4.113c|) . 
fl4.113af . (j4.1i:-ifijl . 04.113/1 ), (|4.115ajl and ()4.1156;i . the terms in the Green's functions 
that involve powers of r or r' do not contribute to the electric or magnetic fields. As 
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we might well have anticipated, only the pure Bessel function terms contribute. (This 
observation was not noted in Ref. |253j .) 

We are now in a position to calculate the pressure on the surface of the sphere from 
the radial-radial component of the stress tensor, 



Trr = - [e{E', + El - El) + ^^{Hl + Hi - H, 



r)\ 5 



(4.120) 



so as a result of the boundary conditions ()4.111|) . the pressure on the cylindrical walls 
are given by the expectation value of the squares of field components just outside the 
cylinder: 



T 



- T 



' ^ ^El + El + -E^ 



r=a+ 
fl — 1 



r=a+ 



Hi + Hi + -Hi 



(4.121) 



These expectation values are given by ()4.9(Jaj) and (j4.96?)jl . where the latter may also be 
written as 

{H{r)H{r')) = — -^{r,r') x V. 



CJ/i 



(4.122) 



It is quite straightforward to compute the vacuum expectation values in terms of the 
coefficients given above. Further details will be supplied in a forthcoming publication. 
The resulting expression for the pressure may then, in a standard manner, be expressed 
after a Euclidean rotation, 



u ^iC, 



A -^ IK, 



so that the Bessel functions are replaced by the modified Bessel functions, 

2 

J^{x')H^{x) -^ —Imiy')Km{y), 
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where 
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A = eyIiiy')KUy) - y'K'Jy)Uy'), 

A = fiyl'^{y')KUy) - y'K'^{y)lm{y')- 



(4.125) 

(4.126a) 

(4.1265) 
(4.126c) 
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This result reduces to the well-known expression for the Casimir pressure when the 
speed of light is the same inside and outside the cylinder, that is, when e^ = 1. Then, 
it is easy to see that the denominator reduces to 

2 = AA = fci-^! [1 - iY[{ImKm)r] , (4.127) 

where ^ = {e — l)/{e + 1). In the numerator introduce polar coordinates, 

y"^ = k^o? + (^0?^ ka = ysm9, Ca = ycos9, (4.128) 

and carry out the trivial integral over 6. The result is 

P = -^^J dyy^ $: -ln(l-e^[y(i^„/.)f), (4.129) 



which is exactly the finite result derived in Ref. |257j . and analyzed in a number of 
papers |259| 1260^ I261J . For ^ = 1 this is the formal result for a perfectly conducting 
cylindrical shell first analyzed in Ref. |25Hj . On the other hand, if ^ is regarded as small, 
and ()4.129|) is expanded in powers of ^^, then the term of order ^^ turns out to vanish, 
for reasons not yet understood j257| 126 if 1^. Recall that the corresponding coefficient 
for a dilute dielectric-diamagnetic sphere ()4.86|) is not zero. 

4-6.2. Bulk Casimir Stress The above expression is incomplete. It contains an 
unobservable "bulk" energy contribution, which the formalism would give if either 
medium, that of the interior with dielectric constant e and permeability fi, or that 
of the exterior with dielectric constant and permeability unity, fills all space. The 
corresponding stresses are computed from the free Green's functions. 



Fl^Kry) = ^G^:^\ry) ^'^ 



e ""' ' ' A'% 



1 / \ I™| 



., I . , . ^. -„.,.... ,-,„,.._, - (4.130) 

2\m\ \ryj 2,1 

It should be noted that such a Green's function does not satisfy the appropriate 
boundary conditions, and therefore we cannot use ()4.12H) . but rather one must compute 
the interior and exterior stresses individually. Because the two scalar Green's functions 
differ only by a factor of /i/e in this case, these are 






A A/7(o) 4- ( y^ - — \ (7(0) 



(4.131) 

while the outside bulk stress is given by the same expression with A' ^ A. When we 
substitute the appropriate interior and exterior Green's functions given in (j4.13Up . and 
perform the Euclidean rotation, u; — > i^ we obtain the following rather simple formula 
for the bulk contribution to the pressure: 

P^ = Tl:f{a-)-Tlf{a+) 

= ^^ E / dW dc [y"Uy')K'M) - {y" + m')iM)KM) 

- yH'^{y)K'^{y) - {y" + m^)I^{y)K.^{y)\ . (4.132) 
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This term must be subtracted from the pressure given in (j4.125|) . Note that this term 
is the direct analog of that found in the case of a dielectric sphere in Ref. |216j - See 
()4.81|1 . Note also that P*" = in the special case e/x = 1. 

In the following, we will be interested in dilute dielectric media, where /i = 1 and 
£ — 1^1. We easily find that when the integrand in ()4 .13211 is expanded in powers of 
£ — 1) the leading terms yield 

oo /•27r -in 



P' 






47r2a4 ^^ Jn ^^7n 27r 

m=—oo ^ " 
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+ -(£ - l)2^^^[/„^(^)K„(y)]' + 0(5-1 
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E / ^yy" 



Uy)KUy) + ^^^j^y[im{y)KUy)r + 0{{e- if) 



(4.133) 
where we have introduced polar coordinates as in ()4.128|) . 

4-6.3. Dilute Dielectric Cylinder We now turn to the case of a dilute dielectric medium 
filling the cylinder, that is, set /i = 1 and consider e — 1 as small. The leading term in 
the pressure, 0[(£: — 1)^], is obtained from ()4.125|) by setting /x = e = 1 everywhere in 
the integrand. The denominator H is then unity, and we get 

^^-^^^ E J^ dyy'Uy)KUy), (4.134) 



m=—oo 



which is exactly what is obtained to leading order from the bulk stress (j4.133p : 

P-P^ = 0[(£-l)2], (4.135) 

which is consistent with the interpretation of the Casimir energy as arising from the 
pairwise interaction of dilutely distributed molecules. In fact, from Ref. |257^ 1262] . we 
know that the van der Waals energy vanishes even in order (e: — 1)^, so we expect the 
same to occur with the Casimir energy, although the latter should diverge in 0[{e — 1)^] 



We now obtain the expression for the 0[(£: — 1)^] term. Because the general 
expression ()4.125p is proportional to £ — 1, we need only expand the integrand to first 
order in this quantity. Let us write it as 

POD OO JY 

dCa / dka ^ — -. (4.136) 



leTr^a^y.^ J_^ ,_.^^_ AA' 



m=— oo 



where we have noted that the (e — 1)^ in S ()4.126aj) can be dropped. Then introducing 
polar coordinates as in ()4.128j) . and expanding the numerator and denominator 
according to 

Ar = iV(o) + (e-l)ivW + ..., AA = 1 + (e - 1)A(i) + . . . , (4.137) 
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the second-order term in the unsubtracted Casimir pressure is given by 

p(2) ^ lizil! r^yy r^Q /^(i) _ ^(l)^(0)^ _ (4_i38) 

Here the correction to the denominator is 

A(i) = yUy)KUy) - V sin^ 9[Uy)KUy)]' + sm^ d K + v') iMKiy) 

-y'sm'eiMKiv), (4.139) 

and the first two term in the numerator expansion are 

iV(o) = -[y^ + m\l - 2sin2^)] lUy)KUy) ' Z/'(l ' ^ ^^^' 0) I'^{y) K'^{y) , (4.140a) 
iV(i) = -i {m' + y^) [y^ + m\l - 2 sin^ 9)] Kl{y)ll{y) 

+ ysin^ ^ [(m^ + y^) + 771^(1 - 2 sin^ 6) - Am^ cos^ O] ll{y)K^{y)K'^{y) 

+ iy2 gin2 ^(1 _ 2 sin^ e)(m^ + y^)ll{y)K'^{y) 

- \y^ sin^ Q {y^ + m2(l - 2 sin^ ^)] Kliy)ll{y) 

+ y^ [sin^ ^ - sin^ ^(1 - 2 sin^ Q)\ I^{y)I'Jy)K^{y)K'Jy) 
+ y=^ [sin^ Q + sin^ 0(1 - 2 sin^ 0)] /;^(y)ir„(yX(y) 

+ \y^ sin^ 0(1 - 2 sin^ 0)/;^ir;^. (4.1406) 

The angular integrals are trivially 

/'27r /'27r o 

/ d0sin2 = 7r, / d0sin^0=-7r, (4.141) 

and then the straightforward reduction of ()4.138|) is 

P^'' = ^^ E rdy{y{y' + m^){2y^ - m^)lliy)KUy) 

m=~OD ^ 

+ 2y\2y' + m2)K^(t/)/„(t/)4(t/) 

- y3(i/^ + m')ll{y)KZiy) - l/=^(2y2 - m')KUy)C{y) 

+ Ay'KUy)K'Jy)Ciy) + 2yH^{y)I'^{y)K'^{y) + y" [I'^{y)K'^{y)f }. (4.142) 
Our challenge now is to evaluate this quantity. 

^.7. Perspective 

I have been working on this problem, on and off, since 1998, when I learned of Romeo's 
proof |262j that the renormahzed van der Waals energy for a dilute dielectric cyhnder 
was zero. Unfortunately, I had labored under a misconception concerning the form of the 
Green's dyadic, which was not in a sufficiently general form until 1 started re-examining 
this problem with my graduate student Ines Cavero-Pelaez this past year. We now have 
a consistent formal result, which only requires some delicate analysis to extract the 
answer. The results, and further details, will follow in a paper to appear later this year. 
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This promises to add another bit of understanding to our knowledge of Casimir forces, 
knowledge that seems to grow only incrementally based on specific calculations, since a 
general understanding is still not at hand. 

5. Casimir Effects for Solitons 

Our discussion throughout this article so far has been confined to idealized boundaries, 
although we alluded to a dynamical basis in the sections referring to the delta-function 
potentials. Of course, from the beginning of the subject, it has always been the goal 
to describe the interactions due to real interfaces, be they constituted of atoms and 
molecules, or due to solutions of the quantum field equations themselves. The most 
natural thing to consider is a solitonic background, where the soliton is a classical field 
configuration which minimizes the energy, and then consider the effect of quantum 
fluctuations around this background field. Perhaps the first physical ideas along this 
direction were presented in the context of the bag model |263t I264( I265| I26b1 lUHl I267t I268J . 
The bag is supposed to represent, semiclassically, the notions of confinement, in which 
within the bag particles carrying color charge (quarks and gluons) are free to move 
subject to perturbative QCD interactions, whereas outside the bag, no colored objects 
can exist. Such a bag picture has never actually been derived from QCD, but it forms 
an enormously fruitful phenomenological framework. Similar pictures can be derived 
from truncated models IM EZOl EZH 12121 EZl EH EZSl • 

Casimir energies have been discussed in connection with the bag model since 1980 
[23 I229| 1^ . (Actually, a zero-point energy parameter was put in the model from 
the beginning.) Unfortunately, no reliable result could be derived because interior 
contributions alone are inherently divergent. Efforts, however, more or less successful, 
were made to extract finite parts |276j , and a summary of some of the phenomenological 
results can be found in Ref. |^. Progress toward understanding the divergences promise 
to lead to more reliable predictions in the near future. 

However, for kink and soliton solutions, reliable Casimir effects have been found 
in a number of cases. The reviews given at QFEXT03 ^4^^ by van Nieuwenhuizen, 
Bordag, and Jaffe are a useful starting point. For example, Refs. |277| I278t I279J show 
that quantum corrections to the mass and central charge of supersymmetric solitons are 
nonvanishing even though zero-point energies of bosons and fermions seem to cancel. 
The Bogomolnyi bound is saturated because there is only one fermionic zero mode. 

Very interesting methods have been presented in the past few years by the group 
led by Jaffe |28U1 I165t 1281] , based on subtraction from the local spectral density (related 
to derivatives of the phase shift) the first few Born terms, which correspond to low- 
order Feynman diagrams, which may be renormalized in the standard manner. For any 
smooth background a finite renormalized vacuum energy is obtained. These methods 
have been used, as noted in section |21 and IH to critically discuss energies and self 
energies of idealized boundaries. In solitonic physics Fahri et al |282j have used these 
methods to compute quantum fluctuations around static classical solitons in Euclidean 
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electroweak theory, which are unstable, in an attempt to find stable quantum solitons. 
(See also Ref. |283j .) No solutions have yet been found, yet some promising nonspherical 
candidates exist. 

Bordag |284j considers the vacuum energy of a fermion in the background of a 
Nielsen-Olesen vortex (string) |285j . The vacuum energy is defined by zeta-function 
regularization, and is expressed in terms of the Jost function, evaluated by using the 
Abel-Plana formula. The quantum correction found in this way, however, is very small. 
It may be that in other cases, such as electroweak strings, the quantum vacuum energy 
might have more physical relevance, even leading to the stability of the string. 

We should also mention that Casimir energies play an important role in lattice 
simulations of QCD. For its role in QCD string formation see, for example, Juge, Kuti, 
and Morningstar |286| I287J and Liischer and Weisz J2l 



6. Dynamical Casimir Effects 

Everything discussed above referred to static configurations. In such a case the concept 
of energy is well-defined, but even then, as we have seen, it is not easy or noncontroversial 
to extract a physically observable effect. When the boundaries are moving, the situation 
is far more difficult. 

In one dimension, the problem seems tractable. We can consider a point undergoing 
harmonic oscillations, and ask what are the consequences for a scalar field which must 
vanish at that point. We expect that the result is the production of real quanta of the 
field. This is the dynamical Casimir effect. However the only reliable results seem to be 
for motions which can be treated perturbatively, or in the opposite extreme, where the 
adiabatic (instantaneous) approximation applies for very rapid changes. 

In three dimensions, the situation is still more challenging. Here we should mention 
the suggestion of Schwinger [213 M, followed up by Eberlein [Ml 1212, Chodos 
ISnSlEnil, Carlson [21312^1, Visser 121111203121321111111112111!, and others, that the 
copious production of light in sonoluminescence |3U2t I236J was due to the dynamical 
Casimir effect, due to the rapid expansion and contraction of a micron-sized air bubbles 
in water. The original estimation that there was sufficient energy available for this 
mechanism was based on a naive use of the cutoff value of | ^ /zu;. An actual calculation 
showed that the energy was insufficient by 10 orders of magnitude |232j . Dynamically, 
photons indeed should be produced by QED by a rapidly oscillating bubble, but to 
produce the requisite number (10^ per fiash) necessitated, if not superluminal velocities 
at least macroscopic collapse time scales of order 10~^^ s, rather than the observed 10~^^ 
s scale jHD]- 

6.1. Fulling/ Unruh/Hawking Radiation 

One regime where definitive results exist for quantum particle production is in the 
general relativistic context. The Moore-Fulling-Davies Effect is the production of 
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photons by a mirror undergoing uniform acceleration |3(J3t l3(J4t l3(J5j . The photon 
spectrum is thermal, with the temperature proportional to the acceleration of the mirror. 
The Unruh effect is very similar |3()6j . If the free equations of quantum field theory are 
examined in the frame of an accelerated observer, with acceleration a, it is found that 
such an observer sees a heat-bath of photons, again with T = ajl-n. (For a precise 
description of these phenomena see the classic book by Birrell and Davies |3U7j .) 

These phenomena naturally are mirrored in gravitational phenomena. The 
celebrated Hawking radiation |3()8j is the production of quanta by a black hole. Energy 
is extracted from the black hole by particle-antiparticle production outside the horizon. 
One particle escapes, while the other falls into the black hole. The resulting thermal 
radiation has a temperature, in accordance with the expectation from the above, 
proportional to the surface gravity of the black hole, or inversely proportional to its 
mass M: 

T = 1.2 X lO^^K ("^y (6.1) 

(Again, see Ref. p7j .) 

Scully and collaborators |3U9j have proposed an experiment to measure the Unruh 
effect by injecting atoms into a microwave or optical cavity, which atoms then undergo 
acceleration. Hu et al jHSl I310J persuasively argue that this experiment will not detect 
the Unruh effect, because the latter does not refer to the radiation produced by the 
accelerated detector (which is nil), Lorentz invariance, crucial to the Unruh effect, is 
broken by the cavity, the thermal distribution of photons in the cavity is not that of the 
Unruh effect, and finally that the injection mechanism will produce cavity excitation so 
that acceleration no longer plays a crucial role. 

For recent work on moving charges, detectors, and mirrors by Hu's group, see 

Ref. ism. 

6.2. Terrestrial Applications 

Most of the calculations of the dynamical Casimir effect have considered scalar fields. 
For example, Crocce et al |312j consider a cavity bifurcated by a semiconducting film 
whose properties can be changed in time by laser pulses, modeled by a time-dependent 
potential 

Ant = -^V{t)6{x)<l)\ (6.2) 

The inhomogeneous wave equation is solved with the time-dependence given as a first- 
order perturbation, with the result that if the film is driven parametrically, that is, 
in resonance with one of the modes of the cavity, particle (photon) production grows 
exponentially. This is in line with the expectations from the 1 + 1 dimensional cavity, 
where if the length undergoes periodic oscillations at a multiple of the fundamental 
frequency ujq = n/L of the cavity 

L{t) = L + AL sin kuot, (6.3) 
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for large times the energy produced is |313| I314| I315| I31(j| ITU] 

Eit) = -'^+^-^^coshkuo^t, A: = 1,2,3,.... (6.4) 

The k = value is the Casimir energy corresponding to ()2.14|) . A numerical simulation 
method for calculating particle production in both cosmological and terrestrial settings 
is given in Ref. |317j . 

One of the few treatments of the 3 + 1 dimensional situation for electromagnetic 
fields is that of Uhlmann et al j318j . who consider a rectangular cavity of length L 
with perfectly conducting walls, with a narrow dielectric slab of width a at one end 
possessing a time-dependent permittivity e{t). The time dependence is still treated 
perturbatively. Only TM modes are effective in producing photons, the number of 
which increase exponentially on resonance, just as in Ref. |312j : 

(iV)(t) = sinh2fS^^ty (6.5) 

where k^ is the square of the transverse wave vector, uj is the resonant frequency, and 
X is the amplitude of the sinusoidal time- varying permittivity, 

—-- = constant + Ysmut. (6.6) 

e{t) 

The challenge will be to devise a practical experiment where this effect can be observed 
in the microwave regime. 

7. Casimir Effect and the Cosmological Constant 

7.1. Cosmological Constant Problem and Recent Observations 

It has been appreciated for many years that there is an apparently fundamental 
conflict between quantum field theory and the smallness of the cosmological constant 
[Slim 13191 I32UJ . This is because the zero-point energy of the quantum fields (including 
gravity) in the universe should give rise to an observable cosmological vacuum energy 
density, 

1 

""cosmo ~ 7^, (7.1) 

where the Planck length is 



Lpi = v^ = 1.6 X 10"^^ cm. (7.2) 

(We use natural units with h = c = 1. The conversion factor is /ic ~ 2 x 10^^^ GeVcm.) 
This means that the cosmic vacuum energy density would be 

Mcosmo~10'''GeVcm-3, (7.3) 

which is 123 orders of magnitude larger than the critical mass density required to close 
the universe: 

Pc = —p^ = 1.05 X 10~^hl GeV cm-3, (7.4) 

SttGat 
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in terms of the diinensionless Hubble constant, h^ = i^o/lOO kms~^Mpc~^ ~ 0.7. 
From relativistic covariance the cosmological vacuum energy density must be the 00 
component of the expectation value of the energy-momentum tensor, which we can 
identify with the cosmological constant: 

{Tn = -ug^" = -^9'"- (7.5) 

[We use the metric with signature (—1, 1, 1, 1).] Of course this is absurd with u given 
by Eq. ()7.3p . which would have caused the universe to expand to zero density long ago. 
For most of the past century, it was the prejudice of theoreticians that the 
cosmological constant was exactly zero, although no one could give a convincing 
argument. In the last few years, however, with the new data gathered on the 
brightness-redshift relation for very distant type la supernovae \6'2l\ I322| 13231 13^ I325J . 
corroborated by observations of the anisotropy in the cosmic microwave background 
|326j . observations of large-scale structure |327j . and of the Sachs- Wolf effect |328j . 
Thus, it seems clear that the cosmological constant is near the critical value, and in fact 
makes up the majority of the energy in the universe, 

^A = A/SttGpc ^ 0.75. (7.6) 

Dark matter makes up most of the rest. Data are consistent with the value for the ratio of 
pressure to energy predicted by the cosmological constant interpretation, w = p/p = — 1 
|329U330J . For reviews of the observational situation, see Ref. |33HI332J . It is very hard 
to understand how the cosmological constant can be nonzero but small. (For a recent 
example of how difficult this problem is to solve, see Dolgov |333j .) 

7.2. Quantum Fluctuations 



In Ref. |334|l335j we have presented a plausible scenario for understanding this puzzle. It 
seems quite clear that vacuum fluctuations in the gravitational and matter fields in flat 
Minkowski space give a zero cosmological constant. On the other hand, since the work 
of Kaluza and Klein |336t 13371 I338J it has been an exciting possibility that there exist 
extra dimensions beyond those of Minkowski space-time. Why do we not experience 
those dimensions? The simplest possibility seems to be that those extra dimensions are 
curled up in a space S of size a, smaller than some observable limit. 

Of course, in recent years, the idea of extra dimensions has become much more 
compelling. Superstring theory requires at least 10 dimensions, six of which must be 
compactified, and the putative M theory, supergravity, is an 11 dimensional theory. 
Perhaps, if only gravity experiences the extra dimensions, they could be of macroscopic 
size. Various scenarios have been suggested |339|l34()|l341j . 

Macroscopic extra dimensions imply deviations from Newton's law at such a scale. 
Five years ago, millimeter scale deviations seemed plausible, and many theorists hoped 
that the higher- dimensional world was on the brink of discovery. Experiments were 
initiated |342t I343J . Recently, the results of definitive Cavendish- type experiments 
have appeared |344| I345| I346| I347J . which indicate no deviation from Newton's law 
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down to 100 ;um. (The experimental constraints on non-Newtonian gravity discussed in 
section ESI are so weak as to be useless in this connection.) 

This poses an extremely serious constraint for model-builders. 

Earlier we had proposed |334j that a very tight constraint indeed emerges if we 
recognize that compact dimensions of size a necessarily possess a quantum vacuum or 
Casimir energy of order u{z) ~ a~^. These can be calculated in simple cases. Appelquist 
and Chodos |3481 I349J found that the Casimir energy for the case of scalar field on a 



circle, S = S^, was 

3C(5) 5.056 X 10-5 

D47r°a* a* 

which needs only to be multiplied by 5 for graviton fluctuations. The general case of 
scalars on S = S^ , N odd, was considered by Candelas and Weinberg j350j . who found 
that the Casimir energy was positive for 3 < A^ < 19, with a maximum at iV = 13 of 
Uc = 1.374 X 10-^/a^. The even dimensional case was much more subtle, because it 
was divergent. Kantowski and Milton (HSl showed that the coefficient of the logarithmic 
divergence was unique, and adopting the Planck length as the natural cutoff, found 

5^iVeven: «^ = ^ln-^, (7.8) 



but ajv was always negative for scalars. In a second paper |351j we extended the analysis 
to vectors, tensors, fermions, and to massive particles, among which cases positive values 
of the (divergent) Casimir energy could be found. In an unsuccessful attempt to find 
stable configurations, the analysis was extended to cases where the internal space was 
the product of spheres |352j . 

It is important to recognize that these Casimir energies correspond to a cosmological 
constant in our 3 + 1 dimensional world, not in the extra compactified dimensions 
or "bulk." They constitute an effective source term in the 4-dimensional Einstein 
equations. Note that because the scale a makes no reference to four- dimensional space, 
the total free energy of the universe (of volume V) arising from this source is F = Vuc, 
so as required for dark energy or a cosmological constant, 

p=-^F = -u, T^'' = -u,g^'', i.e., w = -1. (7.9) 

The goal, of course, in all these investigations was to include graviton 
fluctuations. However, it immediately became apparent that the results were gauge- 
and reparameterization-dependent unless the DeWitt-Vilkovisky formalism was adopted 
|3531 13541 1355| 1356] . This was an extraordinarily difficult task. Among the earlier papers 
in which the unique effective action is given in simple cases we cite Ref. |357j : see 
references therein. Only in 2000 did the general analysis for gravity appear, with results 
for a few special geometries |358j . Cho and Kantowski obtain the unique divergent part 
of the effective action for S = S"^, S'^, and S^, as polynomials in Aa^. (Unfortunately, 
once again, they are unable to find any stable configurations.) 

The results for the coefficient a^- in (j7.8|) are a2 = 1-70 x 10"^, 04 = —0.489, and 
ttg = 5.10, for Aa^ ~ G/a^ <^ 1. Graviton fiuctuations dominate matter fiuctuations. 
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except in the case of a large number of matter fields in a small number of dimensions. 
Of course, it would be very interesting to know the graviton fluctuation results for odd- 
dimensional spaces, but that seems to be a more difficult calculation; it is far easier to 
compute the divergent part, which appears as a heat kernel coefficient, than the finite 
part, which is all there is in odd-dimensional spaces. 

These generic results may be applied to recent popular scenarios. For example, in 
the ADD scheme only gravity propagates in the bulk, while the RS approach has other 
bulk fields in a single extra dimension. 

Let us now perform some simple estimates of the cosmological constant in these 
models. The data require a positive cosmological constant, so we can exclude those 
cases where the Casimir energy is negative. If we use the divergent results for even 
dimensions, merely requiring that this be less than the critical density pc implies the 
inequality (a > 0) 

a> [aln(a/Lpi)]^/^80/im, (7.10) 

where we can approximate (In a/Lpi)^/^ ^2.9. The absence of deviations from Newton's 
law above 100 /im rules out all but one of the gravity cases (S*^) given by Cho and 
Kantowski |358j . For matter fiuctuations only, excluded are A^ > 14 for a single vector 
field and A^ > 6 for a single tensor field. (Fermions always have a negative Casimir 
energy in even dimensions.) Of course, it is possible to achieve cancellations by including 
various matter fields and gravity. In general the Casimir energy is obtained by summing 
over the species of field which propagate in the extra dimensions. 



''^tot 



l^[aan(a/Lpi)+A]^^, (7.11) 



which leads to a lower limit analogous to ()7.10p . Presumably, if exact supersymmetry 
held in the extra dimensions (including supersymmetric boundary conditions), the 
Casimir energy would vanish, but this would seem to be difficult to achieve with large 
extra dimensions (1 mm corresponds to 2 x 10~^ eV). (See, for example, Ref. |333j .) 

That there is a correlation between the currently favored value of the cosmological 
constant and submillimeter-sized extra dimensions has been noted qualitatively before 
|359| I36()| 136 H I362J . A first attempt to calculate the cosmological constant in terms of 
Casimir energies in the context of deconstructed extra dimensions is given by Bauer, 
Lindner, and Seidl |363j . 

In summary, we have proposed the following scenario to explain the predominance 
of dark energy in the universe. 

(i) Quantum fiuctuations of gravity/matter fields in extra dimensions give rise to a 
dark energy, or cosmological constant, oc 1/a^ where a is the size of the extra 
dimensions. 

(ii) The dark energy will be too large unless a > 10 — 300 /im. 

(iii) Laboratory (Cavendish) tests of Newton's law require a < 100 /xm. 
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(iv) Thus, extra dimensions may be on the verge of discovery. If serious hmits on the 
vahdity of Newtonian gravity can be extended down to 10 /im, then we would have 
to conclude that 

(v) Extra dimensions probably do not exist, and dark energy has another origin, for 
example quintessence ^36^ . However, the fact that the rapidly improving data 
favor the cosmological constant interpretation of dark energy |329| 1330) . makes 
alternatives disfavored, since they would generally exhibit time-dependence. 

8. Future Prospects and Perspectives 

In this review we have attempted to present a personal perspective on the progress 
in understanding quantum vacuum energy and its physical implication in the past 
four years. The primary stimulus for the development of the subject has been the 
tremendous progress on the experimental front. This has brought to the fore issues 
that were regarded as arcane, such as the temperature dependence of the Casimir forces 
between metal plates, the meaning of infinities encountered in calculations of quantum 
vacuum energy, and the source of the cosmological dark energy, which it is hard to 
believe does not have something to do with quantum fluctuations, yet is remarkably 
small. At this point, no definitive resolution of any of these issues has been given; 
yet, progress is rapid, and I hope that this status report may help sharpen issues and 
contribute in some small way to the solution of outstanding problems. 

The reader will have noted that this document is far from even-handed. I have 
continued to focus on the use of Green's function techniques as expounded in my earlier 
monograph [SHj. I do not mean to disparage in any way the valuable progress made 
using other techniques, including use of zeta-function methods, Jost functions, worldline 
approaches, and scattering phase-shift formalisms; although I do continue to beheve that 
the Green's function approach is the most physical. I also have focused on topics that 
are of personal interest, so if I have slighted important subjects and researchers, I beg 
forbearance. 

I want to close this review by briefly mentioning a few topics that do not seem to 
fit in the sections above. For example, there has been important work in the subject 
of the Casimir effect in critical systems by Krech and collaborators [36,^ IHTj , in which 
they consider massless excitations caused by critical fluctuations of the order parameter 
of a condensed-matter system about the critical temperature l/Pc- For d transverse 
dimensions, that force is 

/3e^=(t/-l)-, (8.1) 

a 
where A is universal. For an application to thinning of superfluid helium films, see 
Ref. j3()6j . Williams |367j shows that vortex excitations are the source of the critical 
fluctuations that give rise to the critical Casimir force in this situation. 

An acoustic analog of the Casimir effect has been discussed jMHl EEHl EZOI ■ 
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There have been many extremely interesting contributions to cosmological and 
brane-world models. For example, Dowker |371j considers the Casimir effect in nontrivial 
cosmological topologies. A condensate of the metric tensor may stabilize Euclidean 
Einstein gravity in a manner not unrelated to the Casimir effect |372j . And Brevik 
has questioned the meaning of the Cardy-Verlinde formula expressing a bound on the 
entropy |373j . (See also Ref. |374j .) Mazur and Mottola jj3751 I376J have suggested 
that dark energy is quantum vacuum energy due to a causal boundary effect at the 
cosmological horizon - namely, that instead of a black hole, there are three regions due 
to a quantum phase transition (perhaps due to the trace anomaly) in spacetime itself: 
exterior (Schwarzchild) where p = p = 0; interior (de Sitter) where p = —p; and a thin 
boundary shell where p = p. Details of this proposal are still vague; without a detailed 
calculation one cannot tell whether even the vacuum energy will emerge correctly. 

There have been many contributions on Casimir effects in brane-world scenarios, 
for example, Refs. |377| 1378] . 

Finally, we note that everything we have considered in this review has been at the 
one-loop level. Radiative corrections to the Casimir effect have, in fact, been considered 
by several authors. Most of the calculations have been in situations in which there 
is only one significant direction: For QED, see Ref. |379| I38U] . and for A0^ theory, see 
Ref. |38HI382] and references therein. There is an impressive calculation of the radiative 
correction to the Casimir effect with a spherical shell boundary, perfectly conducting 
as far as electromagnetism is concerned, but transparent to electrons by Bordag and 
Lindig P^ . 

So we leave the subject of Casimir phenomena as a work in progress. It is 
clear that quantum fluctuation forces are vitally important both in the very large 
and the very small domains, and that they will play increasingly central roles in 
engineering applications. Thus, the subject is an exciting interdisciplinary topic, with 
both fundamental and technological spinoffs. Thus the uncertainty principle is not just 
about atomic and subatomic physics, but it may control our future, in many senses. 
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